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ABSTRACT 

The design and construction of a position sensitive 
detector, specifically for use in physiological studies, are 
detailed. Such a detector composed of a row of rectangular 
proportional counters, is constructed for studies of the 
movement of Slabelled material along nerve axons. The 
principle of operation of proportional counters, in general, 
is reviewed. 

It is shown that it is possible to optimize operating 
characteristics of a small counter with TeSpeCteco its 
dimensions. Details of a 16-channel prototype multiwire 
Chamber are included, together with the circuitry employed 
in the read-out. 

The application of a multiwire proportional chamber to 
axoplasmic transport studies is discussed. Results ofthis 
are compared to those obtained from similar studies using 
liguid scintillation counting techniques. It is demonstrated 
that the system is readily applicable in thin layer 


chromatography. 
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CHAPTER 1 


INTRODUCTION 


4.1 Introduction 

The development and rapid expansion of nuclear reactors 
in recent years has made available a large number of 
radioactive isotopes which are regularly used in many 
branches of science, either as tracers or as sources of 
radiation. The tracer application is by far the more 
important and is possible because the chemical properties of 
the radioactive isotopes are identical to their stable 
counterparts. One important use of radioactive tracers is as 
radiopharmaceuticals which find wide application in medicine 
and other sciences. In medicine, for instance, they are 
employed to trace metabolic processes, yielding the amount 
of labelled substances present and the rates at which these 
radiolabelled compounds enter and leave the various 
compartments of the body. Another aspect of the isotope 
tracer technique, of particular interest to this work, is 
ifs Bapplicatione ston physiological studies.) Our especitic 
interest is in the study of the movement of material along 


nerve fibers, known as axoplasmic transport. 
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One active area of research in the neurological 
sciences is involved with examining the mechanisms 


associated with the delivery of essential material (e.g. 















+ werTs em / i 
. 7 
nereang he pe. we 
a 
‘apse eats) hah 4 
Lt» sage 7qQ do Gea Ee 
“ol 
co 6 One Din 1 eniasiipg 
“Eee es)! a8 | \7aees ten a) 
; ; ‘a a A076 Aa 
i. rw “yee 
; % oo: thes eee 
: 79 490¢ Gna 
p wh bergen in doabty (cfs ae ee eee 
pa "4 
be ot | oop ers ee aro: te 


ihe TOs ovis sae 
Heept Tonetaieh isi i 


pwyeot tal’ od gdoeing — 


Praia av Ra ou2 ks] Suse 0¥ pygemD 
, 
. . Jedrorgicg, se ,ewpiedves sepaege 
" (atpelbecrten ‘O42 - aide iy | eee avs 7 


‘ A Cu ons i © te aé9 an aot on geed 
, Pap aha y BLAS Wie) Ke wreck «ours. i* 






bal 
7 
_ _ _ a 


is 







proteins) to the axon. It is generally accepted that the 
bulk of protein synthesis takes place in the nerve cell body 
(soma), after which the proteins are transported down the 
axon (see Fig. 1.1). Initially, the study of transport 
phenomenon was performed by damming the nerve, for example, 
by constriction and observing the accumulation of axoplasn 
proximal to the site of constriction! . Definite proof that 
proteins were synthesized in the cell body and then 
transported into the axon was obtained by Droz and LeBlond2. 
Their experiments, performed using {3H]-Leucine, constituted 
one of the earliest applications of a labelled organic 
compound to the study of axoplasmic transport. To date, the 
radioactive isotopes used include 3H, !%*C, 35S, *5Ca and 
32P, The radioisotope tracer technique is now widely applied 
to the detection of both the fast and slow transport systems 
(see below). Detailed descriptions have been given by 
Lubinska3 and Lasek*. 

Measurements of the movement of radioactivity along 
nerves have produced estimates of the velocity of transport. 
For the fast component, which is of specific interest to us, 
the velocity varies from 100 to 2000 mm/day. The slow 
component is limited to a few mm/day. The large variation in 
the fast transport rate may be due in part to the different 
nerve species analyzed and to the methods used to determine 
the transport velocity. 

A typical experimental procedure using radiolabelled 


proteins involves injecting or incubating the nerve cell 
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FIG. 1.1 Simple representation for the fast transport system of 
the site of synthesis of proteins in the nerve soma and 
their transfer into the axon. Protein synthesis occurs 
primarily in the ribosomes of the perikarion. The pro- 
tein produced in the ribosomes is transferred into the 
Golgi complex which loads the protein into the fast 
transport system. The various proteins either travel 
to the nerve endings or are deposited along the axon. 
(After Droz and Koenig, Protein Metabolism of the Nervous 
System 1970. Edited by A. Lajtha) 





bodies with a radiolabelled amino acid PES Gigs OCo a mel hice 
amino acid may be administered either prior to or after 
isolation of the nerve preparation. Following the treatment, 
lasting for a predetermined period of time, the distribution 
of the radioactive protein within the nerve is determined. 
The longer the incubation period, the further the front of 
the distribution moves down the nerve fibers. The rate at 
which the front moves can be determined by elocating the 
front OES the@ activity “for “various “nerve preparations 
incubated over varying time intervals. The distribution of 
the radiolabelled material is normally obtained by liquid 
scintillation counting. Because each preparation yields the 
distribution of the radiolabelled proteins at one particular 
time, several nerves and therefore several animals are 
required. This is in order to obtain a distance-to-time 
relation. The interpretation of the data rests upon two 
major assumptions: 1) that axonal transport is linear in 
time and 2) that all of the preparations are similar, even 


though the nerves originate from different animals. 


3 An alternative technique for the study of axoplasmic 


=. = = soos aoe os = See 


major point of interest arising from their technique is that 
measurements can be obtained from a nerve preparation during 


its incubation using a position sensitive detector, viz. a 
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multiwire proportional chamber (MWPC). 

The nerve preparation is placed in a two compartment 
tray and immersed in a Ringer's solution. In the smaller 
compartment (Fig. 1.2), the spinal cord, roots and ganglia 
(of the sciatic nerve of the amphibian Xenopus laevis ) are 
incubated together with a B-labelled amino acid precursor; 
this compartment is lead shielded. The peripheral nerve, 
after passing through an aperture in the inner wall, is 
accommodated in the larger compartment. Following proper 
placement Ofte hemprepdration= inethes tray, the -aperture) is 
sealed with silicon grease to prevent leakage of the 
radioactive material into the axon (large) compartment. The 
tray is placed on the detector entrance window and the 
transport of labelled proteins monitored for an extended 
period of time via the 8-particles which escape the nerve 
and enter the detector. 

The basic detector components are shown in Fig. 1.2. 
They consist of a collimating system (CS), two drift planes 
(XP and Y®), and an anode plane (Z2P). The collimator, 
provided with four egual voltage steps, confines the 
JOnizing particles (otiginating in (the ,netrve +to (a well 
defined region in space. The electric drift field produced 
by the voltage steps causes the secondary electrons 
liberated between the collimator fins to move towards the 
anode region where each electron gives rise to an avalanche, 
which induces a negative pulse on the corresponding Z-plane 


wire and positive pulses on the X- and Y-plane wires. The 
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Exploded view of multiwire proportional chamber with 
ancillary electronics. CS, collimator system; DL, 
electromagnetic delay line; EW, aluminized mylar entrance 
window; MCA, multichannel analyser; NP, nerve preparation; 
TCA, time-to-amplitude converter; aw anode potential; 

( 


Vpp, drift plane potential; XP; X-(drift) plane; YP, Y- 
(drift) plane; ZP, Z-(anode) plane. (From Snyder et al?®: 
used with author's permission). 
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Charge pulse collected at the anode wire is differentiated 
and the voltage pulse so obtained used to start a 
time-to-amplitude converter (TAC). A summed positive pulse, 
extracted through an electromagnetic delay line (DL) which 
is Capacitively coupled to wires of the X-plane, is used to 
produce a TAC stcp signal. The delay in arrival of the stop 
signal provides information about the avalanche location. 
The system is one dimensional and uses only the X-plane to 
locate the origin of an event along the nerve. The transport 
velocity of the labelled protein along the nerve fibers is 
determined by locating the radioactive front at regular time 
intervals following administration of the radioisotope. 

This new experimental approach to the detection of 
axonal transport has one obvious advantage over the more 
established methods described in an earlier section: it 
yields a distance-time relationship of CLanSPOLIa  LrOmeea 
single nerve preparation. This procedure reduces errors 
arising from changes in ambient condition and, more 
important, iG 2 eliminates errors due to _ physiological 
differences between animals. 

The negative aspects of the detection system just 
described are associated with the special configuration of 
the MWPC required by the axoplasmic transport experiments. 
For example, the need for a collimator, which in the 
experimental set up described above has a centre-to-centre 
spacing of 6 mm, reduces considerably the spatial resolution 


of the detector. An additional disadvantage is the counting 
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loss due to the thickness of the collimator fins. Snyder et 
al.® had to accept a loss of about 25%. The collimator 
dimensions could be modified to reduce these adverse 
effects. There are, however, restrictions on the collimator 
size due to the efficient drifting action obtainable. 

holes that count rate saturation, and therefore efficient 
drifting action, was achieved when the ratio of the hole 
diameter to the voltage plate thickness was in the 
neighbourhood of two. The rectangular collimator used by 
Snyder et al.® had four voltage zones and a total thickness 
of about 16 mm. This would bring the width to voltage zone 
ratio to slightly greater than one, resulting in poor 
collection efficiency of the secondary electrons. 

An additional disadvantage of the MWPC system is its 
background counting rate (15-30 cpm/cm). As reported by the 
authors, the count rate due to the activity present in the 
nerve may be of the same order of magnitude as the 
background, thus necessitating long data collection times 
which degrade the position-time sensitivity of the system. 

The disadvantages just described have been largely 
responsible for the development of a new 8-ray detector for 
axoplasmic transport studies. The detector we propose here 
consists of a series of small proportional counters arranged 
side by side ina linear array. Proportional counters are 
gas filled devices; their principles of operation and some 


general aspects of the counter characteristics will be dealt 
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With in Chapter 2. 

In our particular case each individual counter has a 
rectangular cross section. The dimensions were optimized to 
achieve high counting efficiency to most 8-sources used in 
axonal transport experiments and to minimize the background 
counting rate. These and other aspects of the single counter 
will be presented in Chapter 3. 

The final counter arrangement tested consists of 
sixteen individual counters arranged as a linear array. The 
read-out system employs one amplifier per wire. This 
provides an easy and precise method of identifying the nerve 
region from which ionizing radiation has originated. The 
read-out system and other features particular to the 
combined counter will be discussed in Chapter 4. 

Chapter 5 will contain some experimental results and 
general remarks about the advantages of the detector, and 


also some of its Shortcomings. 
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CHAPTER 2 


PROPORTIONAL COUNTER 


2.1 Gas filled detectors 

The phenomenon of ionization is basic to the detection 
of various types of radiation. When an ionizing particle 
traverses a gas, energy is lost ica the medium due to 
ionization or excitation of the gas atoms or molecules. Of 
interest to us is ionization, which occurs when the energy 
transferred to the medium is sufficient to produce a 
positive ion-electron pair. The energy, W, required to 
create an ion pair for many gases and gaS mixtures is in the 
neighbourhood of 30 eV. Experimental values and methods of 
measurements of the average energy expended by various 
ionizing radiations in a gas medium have been summarized in 
review articles by Valentine and Curran!9%, and Whytet?!. Gas 
filled detectors respond to the ion pairs liberated in their 
active volume by producing an electric signal. Ionization 
chambers measure the number of ion pairs produced per unit 
time, or alternatively the total charge for a fixed period 
of time. Proportional and Geiger devices count the number of 
pulses produced in the counter due to the passage of 
ionizing radiation. 

Gas filled detectors consist, in principle, of two 
electrodes, separated by a potential difference, in a gas 
filled volume. Consider the case of an anode wire lying 


along the axis of a cylindrical cathode which contains the 
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gas mediumt?. At the lower voltages the magnitude of the 
potential difference is an important factor in determining 
the operating characteristic of the counter . The ion pairs 
resulting from the passage of ionizing radiation can 
recombine, in which case the molecules will revert to their 
Original electrically neutral state. However, for 
sufficiently large applied voltages, the positive ions will 
drift towards the cathode where they are reinstated to their 
original neutral state, and the electrons will reach the 
anode causing a current signal to flow in the external 
circuit through the power supply. The dependence of the 
number of electrons produced in the medium on the potential 
difference between the electrodes is demonstrated in 
Figs 2s Vs As the voltage between the electrodes is 
increased, the recombination loss decreases until all the 
pairs produced are collected. At V,/ 4a small additional 
increase of the voltage will + not affect the collection 
efficiency. The device has at this stage reached the 
ionization saturation plateau. This is the region in which 
ionization chambers are Operated.) Tie the Gvoltage is 
increased beyond Vj, the electrons produced by the ionizing 
radiation, under the influence of the electric field, 
acquire enough energy to cause secondary ionization in the 
region near the anode. This process repeats itself until an 
avalanche of electrons reaches the positive electrode. In 
this region, all of the primary electrons will produce an 


avalanche, giving a charge pulse at the anode which is 
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COUNTER VOLTAGE 


Operation characteristic of gas filled detectors. 
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proportional TO the number of avalanches occuring 
simultaneously. For potential differences GLeater  thdns V3, 
thes electrons “acquire enough energy to cause radiative 
collisions with the gas atoms or molecules. This produces a 
continuum of low energy photons which can interact with the 
medium, liberating new avalanches by photo-ionization. The 
number of avalanches produced will exceed the humber of ion 
pairs due to the primary ionizing radiation, and 
proportionality will be lost. at RG a Single avalanche will 
breed more avalanches until the entire gas volume along the 
anode wire becomes ionized. The cloud of positive ions 
counteracts the potential gradient created by the applied 
voltage and ultimately stops the avalanche. 

Beyond V7 follows a plateau which is the region of 
operation of the Geiger-Mueller (G.M.) counter. The pulse 
size obtained in this plateau does not depend any longer on 
the initial ionization, but is a function only of the 
counter characteristics. After the G.M. region (beyond V5) 
the counter enters into a continuous discharge region. 

From the processes just described, it is evident that 
differences exist between the various detectors. Ionization 
Chambers, for example, detect the total number of ion pairs 
produced in the gas volume for any given time. One important 
feature of such detectors is their ability to distinguish 
between different types of ionizing radiation. They are 
frequently used for measurement of radiation dose or dose 


rate. Proportional and Geiger-Mueller detectors, because 

















i fer Va | GAP Sats, Fa ONE er 

| + asp IAF ee : 
; coohn 1 iaiee sites ae wae a 
acne Cagle ire Sjeveeidea™ 


a © 


_ 
Z 


} 7 
| 
499ne enLise GRE PA oe 
wie i2¢ ‘al 1 | a a 


vi = 4 @ ay i aq aA 


cw “SiS 8 2h agate - 
6<€ Sey Se 


| og Soares 


a > - vo v five j } or hee) 


Ld = ay ae et eh: ee 


5 rc Utretedly soe) (begat le 
n 4 Even (rioiuk st2 
b r Arte [eau LA ROR Se Sees 


pr? a an) é is "we ‘as <a toe es GAs 
* ‘s al 4 a ps ( i rel ve ya oe. i , q eet mal 
Prise) “aston. ovenuaw eet (een ett tered gene 


eee AY) Tas or 7 


Soe | waa eae ong yet eee ae 


ioe ae 


a 


14 


they respond to individual charge pulses, are often referred 
to as counting devices. However, since the pulse in the 
proportional region depends upon the number of avalanches 
simultaneously reaching the anode, the proportional counter 
discriminates between the different ionizing radiation, as 


does the ionization chamber. 


2.2.1 Proportional counter: Introduction 

The proportional counter occupies an intermediate 
position between the ionization chamber and the 
Geiger-Mueller counter with respect to principles of 
operation and construction. It operates in the proportional 
region of the gas amplification curve. The counter consists 
of a central electrode (anode) enclosed in a metallic frame 
of any geometry (cathode). The outer case is usually 
provided with an end window for the passage of ionizing 
radiation. 

The principles of operation of the proportional counter 
were first utilized in 1908 by Rutherford and Geiger!3 to 
detect single a-particles. Geiger and Klemperer (1928)!*%, 
using a ball counter, showed that with a suitably adjusted 
potential, the initial ionization produced by each particle 
was magnified linearly by collisions. In addition, they 
found that the amplification factor increased with 
increasing potential up to a critical voltage at which point 
the size of the pulse became independent of the ionizing 


radiation (Geiger Threshold). Barly in its development the 
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proporr ronal Counter iwas “restricted” «to ‘the (detection: of 
highly ionizing particles such as protons and a-particles. 
The lack of electronics capable of resolving charges 
corresponding to a few ion pairs made it impossible to 
utilize proportional counters for the detection of more 
weakly ionizing forms of radiation (e.g. fandy). 

The development of the modern proportional counter 
began in the early 1940's and was directed towards the 
detection of 6 and Vemcadtdt 1On wel nto mMeWwaapplVecat ton 
reguired a complete revision of the old proportional counter 
technology. Certain aspects regarding their operation such 
as gas. mixtures and gas amplification became of increasing 
importance, The earliest complete theoretical and 
experimental investigation of the discharge mechanism in 
proportional counters was reported by Rose and Korff in 
194115, Their results were fundamental towards a better 
understanding of the counter behaviour and the processes 
associated with gas amplification, and stimulated further 
research in the field!6—19, 

In recent years, the increased demand for proportional 
counters for special applications has posed stringent 
requirements on all aspects of the counter technology. Major 
effort has been directed towards exploiting parameters such 
as gas mixture, counter geometry and associated electronics, 
which directly affect counter operations. One innovative 
aspect of the modern proportional counter is the multiwire 


proportional counter (MWPC) 2°. Although multianode counters 
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had been used before, it was customary +o separate their 
active regions by means of one or more shielding wires with 
a conseguent reduction in spatial resolution?!. It was 
believed that for closely spaced wires, their capacitative 
coupling would prevent the localization of an avalanche to a 
Single wire. Charpak et al.29, in their original paper, 
showed that for the detector operating in the proportional 
region, each wire acted as an independent counter. 


The traditional cylindrical counter with a central 


anode is almost exclusively used for studies of gas 
amplification and phenomena associated with avalanche 
formation22-27, The advantage of using a cylindrical 


geometry over other configurations for such work is the well 
understood electrostatic problem. For applications where 
large sensitive areas are required, the MWPC and others 
muitianode configurations have almost entirely replaced the 
cylindrical type. Multiwire proportional chambers find 
application in many areas of intermediate and high energy 
physics for the detection and localization of charged 
particles28-35, As position sensitive detectors, MWPC's are 
well suited to applications in the life sciences for the 
localization of radioactive tracers i biomedical 
studies36-39, Multianode counters are also used in health 
physics for the detection of low energy /Y-rays (<_ 100 
KeV) *0-42, in X-ray astronomy*3, and in radionuclide 


metrology**. 
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22.2 Operation 

The initial process in a proportional counter is the 
ion pair formation in the gas volume. For counter Operation, 
the important step is the subsequent movement of the 
positive ions and electrons towards the oppositely charged 
electrodes. This movement Of Chatges is of  ptimary 
importance to the formation and shape of the pulse. The 
analysis that follows is a précis of the treatment given by 
Wilkinson!2 for a counter of cylindrical geometry, but is 
applicable to counters of other configuration. 

Consider a cylindrical proportional counter (Fig.2.2a) 
with the anode and cathode separated by a potential 
difference. Next, assume the formation of an electron/ion 
pair within the counter volume. The pulse formed at the 
anode results from the collective movement of the positive 
and negative charges in the electric field. Following the 
formation of an ion pair, the potential induced on the anode 
at any time t resuiting from the induced charges “Ga (t) 
and -q_(t), due to the positive ion and electron 
respectively, is 


q,(t) + q_(t) 


P(t) = eee: Gul inne 2-2.1a 


where c is the counter capacitance. At time t=0 the 
potential is zero because the positive and negative charge 
released in near proximity to one another induce on the 


anode signals of the same magnitude but of Opposite 
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Fig. 2.2 Coaxial cylindrical counter: a) cross section b) anode pulse shape. 
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polarity. After collection. of the electron at time t 


iL 
q(t,)=-e and the potential becomes 
de tt.) ae 
Ga 2 gro ele ee ee 
if C 222,15 


At a time t oU%, r>t ), When the positive ion reaches’ the 


cathode, q,(*,) =S0s anda the sporentidlereduces Lo 


ae 2-2.1c 


P(t,) 


In dealing with a large number of ion pairs the problen 
reduces to determining the induced charges -q, (¢) and -q _ (t) 
for the particular counter geometry and types of ions. The 
potential change at the anode results from the contribution 
of all the induced charges. There are some Simplifying 
assumptions which are made in the derivation of the 
expression for the potential for the proportional counter. 
These assumptions are: a) all electrons created in the 
ionization track are collected at the central wire at the 
same time, b) all the positive ion-electron pairs are 
created at the central collecting wire, and c) recombination 
is negligible. The first assumption is reasonable for low 
energy electrons but degrades as the energy increases. The 
second assumption is valid for amplification factors above 
10, Since the avalanche occurs within a few diameters of the 
anode wire. The third assumption is valid for strong 


electric fields. 
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UWergecaictlation, of -qi(t) iand =q: (t): ‘canbe performed 
using Green's reciprocation theorem which states that for an 
isolated system of n conductors, the potentials P due to 
the corresponding charges q; Can be replaced by a second set 
of potentials Ps due to the charges q , to yield the 


equality 


nh 
oy arias >> a BS Ses 
i=1 


We take EF and q! to be the potentials and charges 
respectively under normal operating conditions. Let charge e 
be placed on a very small electrode located between the main 
counter elsctrodes so_ that q, = e; then q and 20 are the 
induced charges to be determined. The corresponding set of 
potentials are taken to be or = PS = 0, leaving 2 
unspecified. The second set of potentials is the one 
existing during normal operating conditions where q, = 0 


and qj and q5 are unspecified. Applying Green's theorem one 


obtains 


q ' a = 
Gia? atad yeh odeP2 Eee:.0 es! 


1?) 
it 
oO 


where CER sd oe + 


follows that 


Hence, solving for eh and q, lke 


Ce ey 

Ne = ee Sie eole 

qd) e Pr = PS) 2-2.4a 
(GENS) 

, = . : 2-2.4b 
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The induced charges qi and q, are determined by the 
potentials under normal operating conditions. In the 
cylindrical counter indicated in Fig. 2.2a, with a central 
electrode of radius 7 and an outer electrode of radius r., 
the potential at the radial position cr. can be obtained from 


the field relation 


eee: VEG 
is in (r,/r,) 


where V is the applied voltage. If we assume that the 


electron is collected at t = 0, from equation 2-2.l1la 


Because of the close proximity of the positive ion to the 
anode at time + = 0, gq, 9) = e and P(0) = 0. Therefore, the 
rise time of the pulse is due primarily to the positive ion 


drifting away from the anode. From Green's theorem 
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2= 27 
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Substituting into equation 2-2.6 gives 


In(r,/r,) 
eis In(r,/r,) 
AGEN WAST Coa GRR 2-2.8 


where r, is the position of the drifting positive ion at any 
time t. 

In proportional counters operated near atmospheric 
pressure the drift velocity of the positive ion, u, iS given 
by 


where K is the positive ion mobility. Thus, from 


ae 
oO 
ay rdr = a; een aR =?) Oh 
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Substituting ee in expression 2-2.8 yields 
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were Vee 
" InG,/t,) /r, ) 





aS a 
AGE eo In; Mites ae j 2-2.11a 
which reduces to 
ny. Ljyaeret ae, 
Ton Gea) 
P(t) = - £ 1 sai 2-2.11b 


Dein (r,/r,) 


At time t2 the positive ion reaches the cathode and the 
potential becomes P(t.) = -(e/c). From equation 2-2.1l1lb one 


can obtain the collection time t>, namely 


ne) 
pee ah y) 2 
t. = Oye a (r, = r,)- 2=2..12 
If we let 
v2 
tec : t 
O 2 ad D 2m, 2 
rey) Tames 
P(t) reduces to 
In(1 + a4) 
1 ae ee meet ne Hs Der eS 


C Z In(r,/r,) 


Pig. 2-2b shows the typical anode ‘pulse obtained in a 


cylindrical counter 
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The number of ion pairs produced in the avalanche, WN, 
is : 
a dul 2-2.14 
where EF is the energy of the primary radiation, W the energy 
required +o produce an ion pair and M the gas amplification 


factor. The expression for the potential thus takes the forn 


(E/W) Me In(l + = ) 
RCC j OO ae 
Ce) é 2 Int,/t,) 222215 
i 
The fast rise time of the pulse obtained in a 


proportional counter is usually exploited by passing the 
Signal through an RC differentiating network, thereby 
reducing the pulse length. To analyze the effect of the 
differentiating network, it is convenient to convert P(t) to 


its normalized form, namely 


Ea 


Z 
1g 
al Z 
i 


where the maximum value that Tt can take is tae ALR 
Knowledge of the differentiated signal allows one to adjust 
the value of the RC time constant so as to yield the maximum 
pulse height while keeping the differentiated pulse to a 


reasonably low width. 
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Fig. 2.3 Differentiating network. 


Given the differentiating network indicated in 


Fig. 2.3, one has that P(t) = q/fc + v(t). Differentiating, 





Gi ne dv (t) 

es C + ir 2-2.17a 
Letting v(t) = if, 

GGoy Buus) dv (tT) a 

a RC ae aie! 2-2.17b 


IORELG.) 265, Ceilssthe capacitance of, the) counter. "9ff the 
counter Signal is decoupled through a capacitor, the 
differentiating capacitance becomes the series sum of the 
fWO. Substituting in z-2.17b for. dP(t )/dt and setting 
1/RkC = S, the first order linear differential equation 


reduces to 


dv (t) 1 1 say 
dt ; 


+ SV(t) =p 
2 In(r4/r,) (t + See, ) 


Deed 
Oia rj/r, 
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The solution for v{t) is found in terms of the exponential 


integral, Fi, and has the form 





2 2 
tr, ry 
Ei US ‘are bd tS = 
S 1 os r3 
V(t) 2? In(r,/r,) 5) 2-2.19 
r 
af 
Exp GS +7) ] 

) 


where 


In the limit of a long differentiating time constant 2c 


(S small), the exponential integral becomes 


Ei (x)) ——_—____ ¥ 4 In x 


me a 10) 


where y is the Buler constant. Then, v(t) reduces to 


1 =, 
V(T) OO in Gl) ere Pt ek 2-2.20 
ss G 2 In(r,/r,) 4 


In particular, the maximum value of v(7) corresponds to the 
maximum value of P(t). 


In the case of very strong differentiation of the 
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counter signal (S very large), 


x 


WSs 1,2! 
Ei(x) = Z [1 + = + 5 + ...] 
x 
and to a first approximation 
1 1 Se 


Pees Deis esT ry ST co a . 2-2.21 
Ty a 

(2 +) (3) 

9 Lo 


BSee5 57 2 7y thes second: term in brackets goes to zero, and 


equation 2-2.21 reduces to 


zs il 1 i Dae Ga) 
Sa ee - 
MASS be Ta eae §. dt ae 
ie 
2 
a 
Since S is positive, the term 
eS. 
a ae 


ec) wD 
S (r,/r5) 


will always be >0 (Equation 2-2.21). Therefore, this term 
will tend to reduce the value .of the term in square 
brackets. The maximum value of V(_) for S>> ri /rj thus occurs 
when the expression 2-2.23 vanishes. The exponential can 


then be expanded in a power series as follows: 
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ae t * al: 


ae ae nee Me ce (-1) . DD ERE 
s(r7/r5) $(r,/r5) =F S 


Applying the ratio test on the infinite series one finds 
that At converges for St =e 1, In) practicey since 
Si(ceeyrdy ) >> 1 the maximum value that V(t) can take is 
Obtained from equation 2-2.22 by letting T = 0 (recall that 
aneeP(T4 5. Thaxt =e) = ri he ). For differentiation such that 
DPA S585 5 Be reequation 2-2.22estill holds, but®an) ethis 


case the time for which V(t) is maximum is determined by1t2 


7 ae 
aaa = In(r,/r,) Te} es 2-2.25 
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In proportional counters, efficient detection of 
ionizing radiation depends largely on the gas amplification. 
In many applications the proportional counter is limited by 
fluctuations in avalanche size due to the probabilistic 
nature of the ionization processes. In general, fluctuation 
in size and shape of the output pulses are acceptable 
wherever the device is used ch) count events. The 


proportional counter, however, is employed in other types of 













fs tad pT vies hy | eas JF 
| value oT << ca 

(wor vos ays or 0 

[ 2 see Ph Vt tee nd an 


aT “9% “doy geet Gee 
L A io A \ —¥e 


ary 4 tity [co 7, a it ae 1 >a 





fej 090 | {hingos Bie Ree eee coeantans lt 


oY AOESOATON. THezolit S443 00ND LSnh se ee pag ot. ja La 


pePe CORAL. WAN Kp Gh? 417 rt nie e bGogeb ot . , Se 


a a | 
aoe 


os, a 
ad 


140w WA 
va a ; =? 
7 





-t 


GO) Petiest ot reranod Lohaksesey sf 
: Ta Pa es PPR 










, 


4 
7 
- | 


20 


measurements such as energy resolution2$°590, particle 
d2-Serimination=* etc.) “in such applications gas gain 
fluctuations must be minimized and accounted for. For 
completeness, a brief analysis of the origin of such 
fluctuations will follow. 

Statistical fluctuations occur in the number of primary 
ion pairs, Nie produced by the ionizing radiation to be 
detected, and in the number of secondary electrons, n=, 
produced in the avalanche initiated by each primary 
electron. The fluctuation in avalanche size becomes’ more 
serious as N decreases, as for example, in the detection of 
soft x-radiation and in ultraviolet radiation. 

The statistics of electron multiplication phenomena © 
have been considered both theoretically and experimentally 
by various authors*5-5!, Curran et al.*© have determined the 
distribution of avalanche size in cylindrical counters by 
releasing Single slow electrons at the cathode 
(photoelectrons). The authors found that the distribution in 
the output pulses due to single primary electrons was best 
fitted by a curve of the form 


L 
t(n) = n* exp(-n) . 2—3.1 


The authors calculated the variance, ar as follows: 
nN 


[e.e) 


if 


oO 


Il 


n P(n) dn f Pin) edn 
O 


tae eo) fan ea 2) w= 
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ao if, na Ptnedn/ fon) dn 
f 


n° =P (7/2)/ T.(3/2) = 15/4 


and 


2 Z -.2 
om =n - (nn) = 3/2 2-3.2 


from which the relative variance is LON Ae) Ramey 5e, lO 
explain the experimental result of Curran etedie 5 SyEne+? 
proposed a model in which the ionization probability is a 


decreasing function of the avalanche size and arrived at a 


Polya or negative binomial distribution of the form 


i 
P(n) = EG ny exp (-n) 2-353 


This particular distribution has a relative variance given 


by 


i 
o 
4+ 
Si ile 
R 
sr 
Bho 
I 
Ww 
= 


for reasonable gas gain in which case n >> 1. The Polya 
distribution reduces to that found by Curran et a 2*6& when. 
b = 2/3 and to the monotonic decreasing distribution, often 


referred to as the Furry distribution, derived by Snyder*s, 
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aaa Onin oN 


oa4 


for b = 1. Byrne concluded that b should take values between 
1 and 0.61, setting a lower limit to the relative variance. 
In addition to the fluctuation in the number of secondary 
electrons produced in the avalanche, there are fluctuations 
in the number of ion pairs My produced by the initial 
jonization. These fluctuations have been analyzed by FanoS2 


who showed that the relative variance should be given by 


2 


oO 
(Ur = aaa 
N N 


where N is the mean number of ion pairs produced and F is 
the Fano factor, which is determined experimentally. 
Alkhazov et al.S3 have measured F values for various gases 
and gas mixtures. 

Provided that there are no interactions between the 
ionization produced by the incident radiation and the 
secondary processes occurring during the avalanche 
formation, probability theory gives the energy resolution of 


the counter as 


where oF is the variance of the measured value of E due to 
fluctuations in the primary DONALZATLONS andeeLln gas 


amplification. The mean number of ion pairs N is related to 
E by WwW, the energy required to produce an ion pair, in the 


following way: 


=m 
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In equation 2-3.7, F and W are characteristic of the gas and 
(On/n) 2is the relative variance in the number of electrons 
in the .avalanche, which can be replaced by the upper and 
lower limits of 1 and 0.61 found by Byrne*7?. 

The gas amplification, and therefore counter Operation, 
depends to a great extent on the gas mixture to be used and 
on the physical characteristics of the counter. According to 
Wilkinson!2, any gas used as an ionizing medium will produce 
gas amplification in a region of high electric field. This 
isapquites true, )showever,.: in a proportional counter, 
requirements such as gain stability, large amplification and 
low working voltage usually limit the choice of gas. Noble 
gases (Ne, Ar, Kr, Xe) are often preferred because they 
provide higher gain. Pure noble gases, however, are very 
seldom used as a gas filling: mixtures of noble and 
polyatomic gases in varying percentages are more frequently 
encountered. In ionizing collisions, noble gas atoms can be 


excited to metastable states of rather long life times (of 
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PNesOCtcmeomylO5*sec) ithe addition lof small “amounts of a 
polyatomic gas favours a faster de-excita+ion of the atoms 
to their ground levels through collisions with the added 
impuritiesi2 94, 

Other phenomena associated with avalanche formation 
which discourage the use of noble gases are processes which 
favour the production of secondary avalanches. These are 
photoelectric emission and positive ion bombardment of the 
cathode?2 26, Both effects are minimized by the addition of 
complex molecules. In fact, the more complex the molecule, 
ther -«dreater the stabilizing effect +.) Charpak et al.256 
suggested the use of gas mixtures of argon with various 


organic quenching agents. Such mixtures are particularly 


useful when the counter is operated in a rather dirty 


environment. 

For well known gas mixtures in common use in 
proportional counters (argon-methane, argon-isobutane, 
argon-carbon dioxide )e gas gains of 106 have been 


reported26 28 , Bemporad et al.SS found that the addition of 
a small percentage of heptane to an argon-propane mixture 
resuited in the suppression of discharge, allowing higher 
operating voltages and consequently higher gain. 

In general, the exponential growth of the avalanche is 
limited by the space charge effect. This occurs when the 
amount of charge liberated in the avalanche becomes. 
comparable to the charge stored in the anode at the same 


locationS®, Bouchier et al.28, adding 0.3-0.6% of Freon 13B1 


on 
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to a standard argon-isobutane mixture (70% Ar - 30% 
isobutane), were able to eee the limit of gas 
amplification set by the space charge phenomenon. The magic 
gas, as it is normally referred to, increases the gas gain 
by about 100, giving an overall gain of 108 without 
entering into the Geiger region. It has been suggested by 
various authors28 57 that the attainment of such large 
amplification is due to the neutralization of the space 
charge effect by the electronegative Freon 13B1. This view 
has been challenged by Lacy et al.S® who claim that rather 
than space charge neutralization, the avalanche beyond a 
certain limit becomes photo-aided without showing the Geiger 
mode behaviour. Gas mixtures containing electronegative 
gases are frequently employed with multiwire proportional 
chambers. These gases reduce the sensitive zone around each 
anode wire to a narrow cylinderS?9?, 

The usefulness of a gas mixture is determined by the 
gas amplification and by the preservation of its chemical 
characteristics over an extended period of time. Gases 
containing hydrocarbons often necessitate regular cleaning 
of the counter because of the polymerization which can take 
place during the development of the avalanche and due to 
occasional sparking. The quality of the magic gas was shown 
to degrade with the age of the mixture, but could be 
restored with the addition of a fourth component, methylal. 
or propylic alcoholsS?. Very often for long term operation, 


the inert mixtures such as argon-carbon dioxide are 
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preferred because they give the counter a longer operating 
life. 

Important from the point of view of gas amplification 
are the requirements placed on the collecting electrode. 


This is so in view of the fact that gas amplification takes 


place near the anode where the field varies as 
tyre le = anode radius). A reduction in the anode size 
allows for a reduction in the working voltage while 


preserving the same field strength. In proportional counters 
the anode radius is usually kept as small as possible (few 
micrometres) in order to obtain high gas amplification. The 
minimum anode size is set by the tensile stregth needed to 
withstand the electrostatic forces6®°®, Gas gain 
non-uniformity resulting from variations in the anode _ size 
along its length has been studied by Dimcovski®!. According 
to the author, the main contribution to gain fluctuation 
arises from local variation of the anode wire and by 
accidental irregularities caused by either impurities 
attached to the anode or sparks. The conclusions were that 
for a fixed anode-cathode distance, the greater the anode 
radius, the higher the relative wire uniformity must be. 

In multiwire proportional counters the anode wire 
spacing s and the electrode gap L play important roles in 
the counter operation. Charpak et al.5*% found that by 
decreasing s, the anode voltage must be increased in order. 
to preserve the same charge per unit length and consequently 


the same amplification. This becomes clear from _ the 
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calculation of the capacitance per unit length c = q/v_ with 
respect to the electrodes. 

For a cylindrical proportional counter, and for that 
matter any Single anode counter, of interest is the effect 
of a reduction in L. Decreasing the gap between anode and 
cathode, the capacitance is increased thereby achieving the 
same gas gain with a somewhat reduced voltage. In 
cylindrical counters the capacitance per unit length is 
given by 1 

eae ei 2-3.8 

#2 i 
For a fixed anode radius r, a reduction in gap length 
implies a reduction of r . The charge stored in the anode 
wire per unit length due to the applied voltage V is given 


by the simple relation 
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CHAPTER 3 


SINGLE COUNTER DESIGN 


Most proportional counters with a single central anode 
have cylindrical geometry. Such devices are preferentially 
used in experimental work concerning phenomena associated 
with their principle of operation. In constructing a4 
detector for axoplasmic transport studies, however, the 
rectangular cross section has definite advantages over the 
cylindrical configuration. Although study of the rectangular 
multianode proportional counter of box-type geometry has 
been reported by other authors*®°, little is known about the 
performance of rectangular proportional counters with a 
Single anode. 

In the construction of a position sensitive detector 
suitable for axoplasmic transport experiments, there were 
three underlying considerations: 

(a) high counting efficiency for 8-radiation, 

(b) narrow field of view of the counter, and 


(Cc) low background counting rate. 


In order to design a detector which meets the criteria 
previously mentioned, we first constructed and tested a. 
number of single wire proportional counters. These devices 


were. rectangular of length 5 cm and various widths and 
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depths. Fach counter consisted of a main body machined from 
a single piece of aluminun, with one side an entrance window 
consisting of 10 micrometres aluminized mylar glued onto the 
counter with conductive adhesive (Eccobond 60-L). The 
detector's internal surfaces were made as smooth as possible 
to reduce possible discharge at sharp points. A 20 
micrometre gold-plated tungsten wire properly tensioned 
(about 15 gms) was placed at the centre of the counter and 


Supported at each end by glass-epoxy blocks. 


HV 


1OM 


COUNTER 
10K 


FIG. 3.1 High voltage network. 


The dielectric fittings in the active region of the counter 
were coated with epoxy resin. One end of the anode wire was 
terminated in the dielectric support and the other end 
connected outside the chamber to the high voltage power 


supply via a 10 M-ohm resistor (see Fig. 3-1). The counter 
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Signal was extracted across a decoupling capacitor and a 
D.C. path to ground provided by a 10 K-ohm resistor. 

To allow for continuous gas flow during operation each 
end of the counter was equipped with a gas line connection. 
A steady flow of approximately 10 cc/min of 93% argon and 7% 
carbon dioxide was maintained by a single stage high 
pressure regulator. The rate was monitored with a Matheson 


flow meter tube no. 600). 


Two types of tests were performed with the various 
counters to examine their detection efficiency for 
8-particles. In the first type, a small source of 14*C was 
placed inside the counter active volume in order that an 
estimate could be made of the absolute efficiency . In the 
second type of test a small source of !*C was placed at 
various positions outside the mylar window in order to 
establish the dependence of the efficiency on position. In 
all studies, the counter pulses were fed via a NB-11 Hamner 
preamplifier into a Tennelec TC-205 linear amplifier. The 
amplified pulses were then passed through a Canberra SCA 
Model 1431 and the output counted with a digital rate meter. 
The SCA lower window was set just above noise level and the 
upper window open to accept the whole sprectrum. All values 
reported here were corrected for background. 

In the first case, small 1*C sources were deposited on 


pieces of mylar which served as the entrance windows of the 
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counters, the sources being centred on the inside of the 
windows. This arrangement yielded, as closely as possible, 
27 counting geometries. Counting rates versus high voltage 
were recorded up to the point of breakdown. To convert these 
rates to absolute efficiencies, the sources were then 
removed, dissolved and liquid scintillation counted. 

The sources were prepared by depositing a drop of 
aceticractd7 sodiumjsalt*(1, 2 ~° 1*c) on mylar “giued . to 
double sided tape. Following evaporation of the solvent, 
sources of about 1 mm in diameter were sectioned out and 
glued separately inside the counters as described above. The 
source thicknesses were estimated to be of the order of a 
few micrometres; this is small enough that self-absorption 
could be neglected. After establishing the HV-count rate 
dependence, the activity of each source was determined by 
liquid scintillation counting. Following removal, each 
source was dissolved in 2 ml of Protosol to which was’ then 
added OF} mi of glacial acetic acid to eliminate 
chemoluminescence . After addition of 10 ml of Econofluor, 
the samples were counted with a Nuclear Chicago Mark I 
liquid scintillation counter. To account for quenching, the 
counting was repeated with 5 microlitres of a standard 
4*C-Toluene added to the sample yielding in such a way the 
counting efficiency. The activity of the source was then 
determined and compared to the count rate obtained with the. 
proportional counter. The error associated with the total 


activity of the various sources was estimated at 10%. This 
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included the error in the standard source (?*C-Toluene) and 
the error resulting from dispensing the 5 microlitres used 
in the internal standard quench correction method. 

Figs. 3.2 and 3.3 show the dependence of detection 
efficiencies on anode voltage for counters with depths of 6, 
12, and 18 mm. Fig. 3.2 was obtained with 4 mm wide counters 
and Fig. 3.3 with counters 3 mm wide. (Measurements were 
also performed using 2 mm wide counters with varying depths, 
however, because of recurring discharge, these counter sizes 
had to be abandoned.) Both graphs show initially a fairly 
linear increase of efficiency with voltage until a plateau 
is reached (with the exception of the 18 mm counters), 
beyond which the counters quickly enter the discharge 
region, as indicated by the rapid increase of the 
efficiency. There are notable differences in plateau size 
depending on the counter width and depth. No plateaus are 
obtained with the 18 mm counters. This can be attributed to 
the recombination of the primary ion pairs due to the weak 
electric field in the region of the counter close to the 
window. (This is better demonstrated by the electric field 
plots for the various counter sizes shown later in the 
chapter.) There is a definite loss in plateau size in going 
from 4 to 3 mm wide counters (from 100 to 60 volts), and a 
decrease in efficiency with increasing counter depth. 
Efficiencies of about 75% were obtained with the 6 mm deep 
counters for either width. In both graphs, the error bars, 


due to counting statistics alone, associated with each data 
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point were maintained at about the 1% level by adjusting 
the counting time, 

In addition to good detection efficiency, it was 
desirable +o construct a counter whose detection uniformity 
across its entrance window would be optimal. Measurements 
were performed employing counters 4 mm wide and depths of 6, 
42, and 18 mm with central anode and placing an external 4*C 
source of half mm in diameter against each counter window. 
The source was produced by screening a wide #*C spot with a 
copper strip supplied with a 0.5 mm circular aperture. The 
activity of the source through the aperture was estimated at 
0.01 microcuries. The three straight lines, a, b and c in 
Fig. 3.4, represent the best detection efficiencies for 
counter depths of 6, 12 and 18 aon respectively with the 
source on the window at the centre of each counter. The 
operating voltage for these measurements was ta yedeed tee 20 
volts, which falls in the middle of the counting plateau. 
For simplicity, we assigned the efficiency of 100% to the 6 
ma deep counter. The graph shows a marked drop in detection 
efficiency as the source is moved from the centre towards 
the edges, and becomes progressively worse with increasing 
counter depth. The error bars in Fig. 3.4 rapresent the 
percentage error which ranges from 1.1% to 1.5% With the 3 
mm wide counters, the drop-off at the two ends should be 
comparable in size because B-particles are confined in a 
narrower region about the anode. 


To determine the usable central region of the counter, 
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FIG. 3.4 Detection efficiency of 
The dashed lines mark the 1% limits. 
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we measured the dependence of the detection efficiency on 
longitudinal position using a counter 3 mm x 6 mm <in 
dimension. In this case the source was moved parallel to the 
anode wire from one end to the other of the counter in steps 
of 5 mm. In Fig. 3.5 the straight line represents the best 
detection efficiency, set conveniently at 100%, while the 
dashed lines enclose the upper and lower limits of 1% error. 
Departure from such limits occurs at 5 mm from each end and 
can be attributed to the loss in field uniformity due to 
counter end-effects. 

According to these preliminary tests, the 3 mm x 6 mo 
configuration was found suitable for the construction of the 
multiple counter to be described in Chapter 4. However, the 
lack in counting plateau size tends to discourage the use of 
such counters because the detection efficiency would 
unfavourably depend on voltage fluctuation with time. 
Nevertheless, given a stable power supply, the change in 
efficiency can be confined within acceptable limits (less 
than 0.1%); this is true in our case where a FLUKE Model 
408B is used which has a stability of 0.02% per day after 
warmup. A compromise between plateau size and Spatial 
resolution has to be accepted in axonal transport 
measurements. This is achieved by using a smaller counter 
width, thereby increasing Spatial resolution and 
compensating, to a certain extent, for the smaller counting 
plateau. 


With the multiple detector the depth of the individual 
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counter was increased to 20 mm to ease construction 
problems. The optimal 3 mm x 6 mm dimensions were preserved 
by displacing the anode wire off centre to 3 mm from the 
detector window. The electrostatic study of such a counter 
geometry, performed in section 3.3, indicates no major 
change in field configuration suggesting a close similarity 


between the 6 mm and 20 mm deep counters. 


3:3 Electrostatics of the counter 

Proportional counters of rectangular cross section are 
preferred over cylindrical ones when an external source must 
be viewed using a side window configuration; in addition, 
large sensitive volumes are easily achieved with a single or 
with an array of such rectangular counters. Inherent 
disadvantages are that the field distribution is not 
isotropic and that equipotential lines become increasingly 
distorted with increasing distance from the anode, therefore 
affecting the avalanche formation. 

The electrostatic problem in two dimensions for the 
rectangular counter was dealt with in detail by Tomitani®2. 
The potential is written in analytical form using a 
conformal transformation of the Jacobian elliptic function. 
The method, however, is particular +o an anode _ placed 
symmetrically in the chamber. Urbanczyk and Walisorski®3 
approached the problem using the principle of finite 
difference approximation. This technique is a lengthy one, 


and in addition, the error in the approximation is affected 























dohteterenos esep 4S te Sees er 

a a onotoneet} tht hi ae | 
“iy Ghee) eae ricdnihinae : 

‘abt an? «atbals spaniel i 

oY ‘haeveisay Tadapoke 

(ipetso> bert iat «pee a4 

nt ‘ah eo } of) coaeeed ; : 


~ 


wj6bus iY Fey og lesa ee 
1) 2TAr ayes (2RGh Sn QeRe : ‘ 
ss Je6e ta en bi oa 
nev obj ® Babes sueeev ede -_ 
jxeei4 \ be ee gaa a tDibigd’ F BA ‘(joes sptet . 
ieprerwa bile 22. teeta, oS (Tee, ' 
77 a4 : noone vires : 
biuhbi-<s (76 tenth ailt tia att toes - 
‘ wla kote) rer ieraafalen 
noeeed dda rs a wn 
ig =a* ss ruil ee ee ee ie 
PS Ped tiaeT oo dociete bons wal > ho) haw sec néoe . 
. : ) ae 
ae ee 
m4 Poe 


b 7 7 y an : s 


a yoa- e507... Levatyiegn owe 


- 


48 


by various factors. 

For the rectangular counter considered here it was 
found more convenient to use the technique based on the 
integral representation of potential. This method was used 
by Symm®S to find equipotential contours for semicircular 
cross section and pill-box type detectors. In this 
representation there are no restrictions on the shape of the 
cathode and location of the anode. The numerical solution to 
the system of linear integral equations is easily obtained 
and contour plotted using a computer. Lines of equal field 
strength are obtained by numerical differentiation of the 
potential function. Details of this technique are contained 
in Appendix A. 

This numerical method was applied to a square counter 
with an anode of 20 micrometres diameter placed at the 
centre. Both equipotential lines and lines of equal field 
strength were obtained and plotted (Figs. 3.6a and 3.6b). 
The equipotential lines (Fig. 3.6a) are plotted for 
potentials from 5% to 60% of the anode, in increments of 5%; 
the potential decreases with increased distance from the 
anode. It is obvious from the plot that there is some 
Similarity with the circular case. The distortion near the 
corners is obviously due to the geometry of the counter. 
This corner effect shows up more clearly in Fig. 3.6b. The 
field strength drops more rapidly from anode to corner than 
it does from anode to side of counter. The field strength 


difference between the closely spaced lines in the corner is 
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(b) 





FIG. 3.6 Square counter with central anode 
a) equipotential lines; b) magnitude field lines. 
Refer to text for additional details. 


(a) 





FIG, 3.7 Rectangular counter with central anode and 
2:1 length:width ratio a) equipotential lines; 
b) magnitude field lines. 
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10% of that between the wider spaced lines closer to the 
anode. 

A Similar analysis was applied to a rectangular 
counter. Once again, the anode is located symmetrically in 
the chamber. The results for a counter whose height to width 
Ea .Omete wee are sillustrated in Figs. 3./4a and 3.76. The 
equipotential again exibits a smooth decline with distance 
from the anode. However, as in Fig. 3.6b, Fig. 3.7b shows 
clearly the anisotropy of the field strength. Results for a 
counter whose ratio of dimensions is 4:1 are illustrated in 
Figs. 3.8a and 3.8b. A disconcerting aspect of the 
anisotropy of the field is the weakness of the field at the 
corners. Ton pairs produced in these regions of the counter 
have relatively long drift times, therefore enhancing the 
probability of recombination. Another important 
consideration to be made is that electrons liberated in 
these regions are more likely to be captured by 
electronegative impurities present in the detector than if 
they were released in regions of higher field intensity. 

To preserve the field strength towards the upper end of 
the counter which holds the entrance window, the anode was 
displaced off centre. This particular configuration was 
necessary for the construction of the multiple detector and 
described in the following chapter. Figs. 3.9a and 3.9b are 
the result for a counter with sides of 2:1 ratio and anode 
displaced 3/4 of the height. In the upper quarter of the 


counter the egquipotential and absolute field resemble the 
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(a) 














FIG. 3.8 Rectangular counter with central anode and 
4:1 length to width ratio a) equipotential lines; 
b) magnitude field lines. 





ae 


(a) (b) 
LY) 





FIG. 3.9 Rectangular counter with off-set anode. 
a) equipotential lines; b) magnitude field lines. 


3 


Square counter with centrally located anode. The lower part 
is similar to the rectangular counter (See Fig. 3. 29a /and 
Fig. 3.8a). This almost field free region can be extended 
without significantly affecting the upper active region of 
the counter as is shown in Figs. 3.10a and 3.10b where the 
Sides are in 4:1 ratio and the anode located at 3/4. of. tts 


height. 


The electrostatic study of a counter with off-set anode 
Showed that the equipotential contours towards the window 
were little affected as the counter depth was increased, 
provided the anode location remained unchanged. It was 
reasonable to assume that the detection characteristics of 
the 3 x 6 mm2@ counter with central anode could be preserved 
if the height was increased to 20 mm, keeping the other 
parameters unaltered. Using a counter of dimensions 3 mm by 
20 mm (width-depth respectively), with the anode 3 mm _ from 
the window, we determined the detection efficiency as 
described in section 3.2. Fig. 3.11 shows the efficiency of 
B-particle detection versus anode voltage. The size of the 
counting plateau remains at about 60 volts with an 
efficiency ranging from 65% to 72% between the lower and 
upper ends of the plateau. 

The relative detection efficiency across the counter 
window was again determined adopting the procedure used with 


the different 4 mm wide counters (See section 3.2). In 
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FIG. 3.10 Rectangular counter with off-set anode and 
4:1 lenth:width ratio. a) equipotential 
lines; b) absolute field line. 
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Fig. 3.12 the straight line represents the optimal detection 
efficiency, set once more for SiNplicatyee dt.—6100%,- “and 
corresponding to anode voltage of 1150 v. The Magnitude of 
the error bars accounts for the fractional error, eHich 
Tanges trom 1.6 to 1.8%. 

A typical 8-spectrum of an external 14C source obtained 
with the 3 mm by 20 mm counter is shown in Pig~ 36 13. THe 
amplified pulses were, in this case, fed into an ORTEC 6240 
NCA for pulse height analysis. The spectrum shown was 
obtained at the plateau voltage of 1150 Vv. The system was 
calibrated using the 5.9 Kev Mn K x-ray of S5Fe and assuming 
a linear relation between channel number and energy. The 
Spectrum itself represents the response of the counter to 


the distribution in energy deposition within the gas volume. 


In axoplasmic transport studies, the activity to be 
measured depends on the uptake of labelled amino-acid 
precursor by the nerve soma and the amount of active 
material which is actually transported down the nerve. 
Indications are that only 50-200 ppm of the original 
activity is transported into the axons. Due to the geometry 
of the experiment, only a small fraction of this activity 
(0.2 - 2% ) will actually be detected by the counter. As an 
example, for nerves incubated with 0.5 mCi of 1%C or 355 in 
a2.5 ml of Ringer's solution, the expected counting rate 


will be of the order of a_ few counts per second. This 
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FIG. 3.12 Detection efficiency of 1*C across the counter window 
(3 mm by 20 mm). 
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necessitates a low background counting rate. 

The background arising from isotopes of the counter 
material and gas, such as 26Al and !*C, are negligible . 
Background due to cosmic radiation should be negligible 
mainly because of the small counter volume; measurements of 
background performed with and without lead shielding showed 
a difference of 4%, attributable in part to cosmic 
radiation. Fig. 3.14 shows the background spectrum of a 3 mm 
wide and 20 mm deep counter ( without lead shielding) taken 
over a period of 4000 sec. The electronic noise pick-up is 
indicated by the sharp rise at the lower end of the spectrum 
and is. normally reduced by raising the lower level 
discriminator. Threshold adjustment is needed whenever the 
noise level is of the same order of magnitude as or greater 
than the actual activity under measurement . In our case the 
threshold level was kept as low as possible to minimize the 
loss of §-spectrum at the lower energy end. Operating under 
typical conditions, a background counting rate of 6-8 cpm 


for a 3x20 mm2@ counter was observed. 


One important aspect of proportional counter operation 
is the signal amplification resulting from the ionization of 
the gas mixture by the primary radiation and the following 
formation of the avalanche at the collecting electrode. The 
charge constituting the avalanche depends in general on the 


gas mixture and on the physical characteristics of the 
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FIG. 3.14 Background spectrum for a 3 mm by 20 mm 
counter. 
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counter. With reference to our particular counters, of 
interest is the effect that the anode-cathode spacing has on 
gas amplification. Mixtures of argon and carbon dioxide are 
often used as the interacting gas medium in proportional 
counters. The gas adopted in our work had a composition of 
93% argon and 7% carbon dioxide, and was obtained from 
Matheson Co. Inc. 

The counters used for gas gain measurements had 
rectangular cross sections with depths of 6, 12 and 18 mn, 
and widths of 3 and 4 mm. The anode in all cases was located 
at the centre of the counter. Gas flow was maintained 
constant throughout a measurement using a Single stage high 
pressure regulator. The flow rate was monitored with a 
Matheson tube flow meter. Each detector was first purged by 
a fast flow of gas and then allowed to stabilize at a low 
rate. The equilibrium condition was checked using a 
collimated S5Fe source by noting the shift in peak position, 
using a multichannel analyser ( ORTEC 6240A). The anode 
potential was provided by a FLUKE high voltage power supply 
(Model 408B) with an accuracy of 0.25%. 

The counter pulses were fed through a charge sensitive 
preamplifier (ORTEC 109A) into a Tennelec TC-205 linear 
amplifier having a shaping constant of 3 microseconds. The 
Signal was then analyzed using the MCA. To calibrate the 
system, voltage pulses with a rise time of 20 nanoseconds 
and a decay time of 200 microseconds were fed through a 


known capacitor into the amplifying system and the pulse 
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height analysed by the MCA. A series of calibration curves 
showing the relationship of charge to channel Number for 
various amplifier gains was obtained. 

A collimated SSFe source emitting 5.9 keV Mn K x-rays 
was placed against the counter window. The Charge resulting 
from the absorption of the photons at different anode 
voltages was fed through the amplifying system. The location 
of the photo-peak permitted the determination of the total 
charge produced by the Ny ion pairs where Rep (5900 eV) MW. 
For our particular gas mixture the average ionization 


energy, W, (26.6 eV) was computed using the approximate 


relation suggested by Meyers¢6; 


In the above expression P is the relative partial pressure 
LOE the i-th gas component and Wy the corresponding 
ionization energy for the pure gas, obtained from Whyte!!, 
This technique cf comparing a known amount of Charge to 
the total charge generated in the detector has been used by 
other investigators!9»66 69, Its inadequacies stem from the 
fact that the rise time of the charge pulses from the 
precision pulser-capacitor, and the detector may be 
Significantly different, in which case the pulse-shaping 
networks in the preamplifier and main amplifier will respond 


differently. Thus, two identical quantities of charge 


derived from the two devices will produce pulses which are 
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different at the Output of the main amplifier. In absolute 
gas gain measurements the attenuation produced by the 
Shaping networks must be taken into account. Otherwise, the 
calibration procedure is quite acceptable in cases such as 
ours where aproximate gains are required, as in 
demonstrating the dependence of gas gain on counter size. 
Fig. 3.15 is the graphical representation of the number 
of electrons produced by an ion pair versus anode voltages 
for a number of counter. sizes. In the voltage ranges 
considered the results obtained compare favourably with 
values reported in the literaturesé, The gas gain-voltage 
characteristics for the three counters of 3 x 6 mm@, 3 x 12 
mm@ and 3 x 18 mm2 (width x depth) are presented in 
Fig. 3.15a. The gas amplification for the first two of these 
are almost identical to one another, as indicated by the 
overlap of the two lines. The gain for the 18 mm counter is, 
however, consistently lower. Similar results were obtained 
with 4 mm wide counters and are illustrated in ELgee oe toe 
In this case, there is again a consistent decrease in gain 
with increasing counter depth. The dependence of gain on 
counter width is demonstrated in Fig. 3.16. Line (a) shows 
the gain in a 3 x 6 mm2 counter while line (b) is obtained 
for a 4 x 6 mm2 counter. The drop in gain with increased 
width for fixed anode voltage may be attributed to the 
reduction in capacitance per unit length in the wider 


counter. 
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WIDTHS: (6) 3 mm 


(@) 4 mm 





800 900 1000 1100 1200 1300 


ANODE VOLTAGE 


FIG. 3.16 Gas amplification for 6 mm deep counters. 
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3:/ Discussion 

The objective of the study described in the present 
Chapter was to test various counter sizes of rectangular 
cross section and to establish the optimal dimensions with 
regard to detection efficiency and background. Such a 
counter would represent one unit in a larger eee 
sensitive detector applied to axoplasmic transport 
measurements. The detection characteristics of the 
rectangular counter were compared to the position sensitive 
detector, a MWPC, employed in the axoplasmic transport 
measurements of Snyder et al.&. Adopting a counter size of 3 
x 20 mm@ for the multiple detector (see Chapter 4), the 
spatial resolution was improved by a factor of 2 as compared 
to the previous system. The average detection efficiency, 
that is, the average efficiency at any point of the active 
counter region including the collimator, was increased from 
48% to 60% with the new detector size. The results of the 
gas gain measurements indicated that the counters were in 
fact operating in the proportional region, and that no 
Significant increase in gas gain could be anticipated prior 
to entering into the Geiger or discharge regions. In 
addition, the new counter had the advantage of a slightly 


lower background. 
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CHAPTER 4 


MULTIPLE DETECTOR 


A multiple detector was constructed based on the 
results obtained from the single counter studies. Some 
Single counter features were modified to overcome multiple 
counter construction difficulties. Improvements were 
introduced, with particular regard to discharge phenomena. 
The multiple counter is composed of a linear array of 16 
Single chambers separated from one another by thin walls, 
referred to as "fins". A major effort was directed to the 
assembly of the various counter parts. In particular, the 
counter fins and wire plane were mounted in a manner which 
allowed for ease of replacement. For the continuous gas 
flow, we provided the counter with two manifolds (inlet and 
outlet) while allowing for separate gas apertures to each 
counter. 

After considering various readout systems, the "one 
amplifier per wire" approach was chosen. The amplifier used 
is similar to that described by Tarlé et al.69 , with some 
changes introduced, mainly in the logic circuitry. It yields 
a non diinear amplification with output of adjustable width 
and fixed height. Boards which contain 16 such amplifiers 
and associated logic circuitry were made available to us by 
the TRIUMF, groupyoft the NRCG atuthe Us of “A. 


After amplification, each counter signal is assigned a 
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binary code so that data from each counter can be stored in 
a fixed location of a multichannel analyzer (MCA) memory. 
The data stored in the MCA is regularly dumped onto paper 
tape for subsequent analysis. 

The multiple counter of 16 channels plus amplifiers and 
encoding system represents a basic unit which can be easily 


expanded. 


4.2 Counter con 
The multiple detector consists of 16 compartments 
disposed in a linear array. The individual chambers are 3 x 
20 x 5O mm3 (width x depth x length) and are arranged side 
by side. The compartments are separated from one another by 
thin metal walls. Each chamber is a proportional counter 
with an anode wire located at width median and at a distance 
of 3 mm from the counter window. The multiple counter has an 
outer aluminum structure which provides a gas tight 
compartment, and in addition, electrostatic shielding. 
Inside the main metallic frame is a system of four 
planes which hold the fins and maintain the wires in the 
same plane. Fig. 4.1 shows the four planes which are made of 
Nema G-10 glass epoxy board. Plane (a) has centre spaced 
vertical grooves as shown. They provide spacing and support 
for the fins, which are inserted from the top. The top 
horizontal channels surround the anode wires, which run 
along the centre of the counters. The function of the 


horizontal channels is to increase the dielectric surface 
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path length to ground, thus reducing surface leakage. The 16 
20 micrometres gold plated tungsten anode wires are mounted 
under slight tension (about 15 gms) on plane (b). Plane (b) 
and the spacing plane (c) fit around the elevated parts of 
piane (4) aS illustrated in Fig. 4.1. Plane {d), with the 
upper face copper plated and kept at ground potential, 
provides the final spacing in the counter and helps to 
Maintain the window at ground potential by simple contact. 
Those non metallic surfaces enclosed in the active regions 
of the counter which required machining were coated with 
Eccobond 24 epoxy resin. This procedure helps to maintain 
clean and smooth surfaces, preventing leakage of gaseous 
impurities into the counter and reducing discharge at sharp 
points. 

The four planes are fixed to a metallic frame, shown as 
section (a) in Fig. 4.2. The elevated part of section (a) 
has equally spaced slits for the counter fins, and in 
addition, small apertures at each end for continuous gas 
flow. These apertures are connected to the main ducts 
indicated by the tube fittings on each side of the counter 
base. 

The mid section in Fig. 4.2 is provided on each side 
with suitable channels for the housing of O-rings. The front 
end contains a Vero edge connector (part no. 2231-1) which 
provides the contact between the wire plane and the high 
voltage power supply. The connector, made of Diallyl 


Phthalate, has a capacitance between adjacent contacts of 1 
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FIG. 4.2 Multiple detector: 


metallic frame. 
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pF; the proof voltage D.C. at sea level is given at 1,350V. 
The connector is attached to a high voltage network located 
in an isolated compartment of the unit housing the amplifier 
board and the encoding circuitry. Each anode wire is coupled 
through a 10 M-ohm resistor to a common power supply, as 
Shown in Fig. 4.3. The 16 outputs, decoupled through 0.05 
microfarad capacitors, C , are connected to the inputs of 


their respective amplifiers. 





OK 





FIG. 4.3 High voltage network for the multiple detector. 


The top section, (c), has a rectangular aperture for 
the 20 micrometres aluminized mylar counter window. The 
window is glued to the frame with Eccobond 60-L conductive 
adhesive to provide ground connection. 

When the counter is assembled, the three sections are 


bolted together to achieve a gas tight compartment. The 
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metallic case is kept at ground potential and provides 


electrostatic shielding for the counter anode wires. 


—eeees 90 Sree ae ee a es ee re ee ee Se ee es 


The main function of the fins is to delimit each 
counter volume and to prevent the ionizing radiation from 
penetrating more than one counter. They are metal strips of 
predetermined thickness, grounded by simple contact to the 
metallic case at the base of the counter which holds them in 
place. The fin thickness is determined by the energy of the 
radioactive source used in axoplasmic transport experiments, 
in our particular case, 8-emitters whose maximum energy, En! 
varies from 0.154 MeV for 1%*C to 1.71 MeV for 32p, 
Nevertheless, their thickness must be kept as small as 
possible to prevent excessive screening of the nerve from 
the detector. 

The absolute range of beta particles in a given 
absorber is not a well defined quantity. The experimental 
determination of range is not easily performed because the 
measurements are made at practically zero intensities. To 
estimate the maximum range, Rue of some common g-sources, 


we used the empirical energy range relation obtained by Katz 


and Penfold79, 
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where n = 1.265 - 0.0954 1n En. The maximum energy, E,, is 
expressed in MeV, and the range is independent of the ‘type 
of absorber. The above equation is valid for both 
monoenergetic electrons and continuous f-emitters, and holds 
true for energies between 0.01 - 2.5 MeV. 

Table 4.1 gives the thickness in cm x 1073 for various 
metal absorbers and different $-sources. The maximum range 
is obtained from equation 4-2.1 using the known £, values 
for four common ®@sources used in axoplasmic transport 
studies. It is obvious from the table that for sources with 
ranges less than or equal to that of *5Ca, an Al fin of 
thickness 0.5mm would provide effective shielding between 
neighbouring wires. In addition, since the counter was 
devised with easily removable fins, it is possible to use 


the same detectors with a 32P source merely by replacing the 


aluminum fins with tungsten fins of the same thickness. 


One important aspect of the multiple counter which must 
be dealt with in detail is the gas flow. In detectors of the 
continuous flow type, the gas mixture must move uniformly 
throughout the counter. This allows for a proper washout of 
electronegative impurities such as O, , and maintains a 
uniform detection efficiency throughout the counter's active 
volume. The rate of flow of the gas mixture in each chamber 
should be such that even with a slight change in the flow 


rate, the detection effeciciency should remain unaltered. 
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The counter consists mainly of 16 compartments, each 
provided with an aperture at each end connected to the main 
gas lines. One line leads to a pressurized gas cylinder, 
containing the gas mixture of 93% argon and 7% carbon 
dioxide, via a single stage high pressure regulator anda 
flow meter. The second line is open to the atmosphere 
through a long narrow exhaust tube, which keeps the counter 
at a slightly elevated pressure, and causes a slight bulging 
of the thin counter window. The pressure gradient between 
the two ends of the long tube maintains an outflow of gas 
which in turn prevents an inflow of gas impurities. 

The optimum flow rate of the gas through the multiple 
detector was determined by measuring the counting rate with 
respect to the total gas flow. The supporting electronics 
consisted of a Hamner NB-11 preamplifier and a Tennelec 
TC-205 amplifier with a total gain of about 1000. The 
amplified signal was fed into an ORTEC digital ratemeter 
(Model 434) through a Camberra SCA (Model 1431) with the 
upper window set to cover the whole spectrum. This 
arrangement of the recording system was used throughout in 
this section. Prior to any measurement, the detector was 
outgassed for several HOULS! “at dee rates Of 10s stc/min. 
Following each change in the flow rate, enough time was 
allowed for the flow to stabilize. Fig. 4.4 was obtained by 
placing a small 1*C source (about 1 mm in diameter) against 
the counter window, centered with respect to one particular 


Channel. Each data point was obtained from a set of 5 
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independent measurements. The 100% level was set to 
correspond to the maximum mean count rate obtained from the 
group of measurements. The error bar associated with each 
data point gives the percentage standard deviation. Such 
deviations range from +0.3% to +0.7%. 

Fixing the flow rate at the value of 80 cc/min, the 14c 
source was moved along the axis of one single counter, 
thereby determining the detection uniformity in the 
longitudinal direction. Fig. 4.5 shows the dependence of the 
count rate on the position of the source with respect to the 
two gas apertures. The mean of the set of measurements was 
set at 100%. The dashed lines above and below the mean value 
mark 1% deviation from the mean. 

The response of each of the 16 counters in the multiple 
detector to the 8-source is shown in Fig. 4.6. Once again, 
the total gas flow rate was kept at 80 cc/min. The source 
was placed in contact with the window at the centre of each 
counter. Additional measurements were made at 5 mm from each 
end of the 16 counters. It was found that changes in count 
rate were within +1% of the mean. AS in Fig. 4.5, the 
dashed lines in Fig. 4.6 represent 1% deviation about the 
mean. 

From the results shown in Fig. 4.5 and 4.6, one can 
conclude that the variation in detection due to gas flow is 


within an acceptable limit. 
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The; most critical aspect of the read-out system is 
related to the low detection threshold. In detecting weakly 
ionizing radiation, such as B-particles, this threshold may 
correspond to no more than a few ion pairs. With good gas 
amplification, this would amount to only a few hundreds of 
microvolts. 

Following the development of MWPC's Charpak et al.71 
described methods which could lead to Simple and inexpensive 
read-out systems. Two of the techniques described by the 
authors are noteworthy because of the Simplicity of the 
principle of Operation which leads the way to more 
sophisticated methods: the lumped delay-line technique and 
the pulse division method. In the first of these, the wires 
are separated by inductances. The time difference between 
the arrival of the pulse at the Opposite ends determines the 
wire carrying the signal. The pulse division method uses a 
purely resistive separation to produce a detectable pulse 
attenuation while preserving the pulse Shape. To remove the 
effect of the spread in the initial pulse heights, the 
location of the wire is determined by the ratio of the 
pulses appearing at the opposing ends. 

More recently, a method of reading the position of an 
ionizing event making use of commercially available delay 
lines was developed by Rindi et al.72 and Grove et al.73 The 
principle of operation is very much similar to the lumped 


delay line technique suggested by Charpak. The construction 
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and characteristics of an improved version of a delay-line 
has been given by Grove et alwi*” oP ine generals) »the 
positioning accuracy which can be achieved with these delay 
lines depends on the amplitide of the pulse above electronic 
noise and the polarity and Spacing of the wire plane which 
is read-out. For the electron collecting plane the accuracy 
amounts to half the anode wire spacing. A negative aspect 
which makes the use of delay lines with our detecting system 
less attractive is the attenuation produced by the 
capacitative coupling between the wire plane and the delay 
line. For good delay lines a coupling efficiency of 18% has 
been reported?’s, 

In our read-out system we adopted the one amplifier per 
wire method; the amplifier, trigger and memory are described 
by Tarlé et al.©%. One important reason for choosing this 
particular circuitry is the low voltage threshold of the 
amplifier and trigger, normally .5mV as reported by the 
authors. One practical aspect of the amplifier, trigger and 
memory circuit, in addition to the low cost, is that it 
amplifies the counter signal to a level which permits the 
handling and interpreting of the information with minimal 


additional electronics, 


The amplifier used with our read-out system was. 
originally designed by Tarlé et al. . It utilizes a Motorola 


NC 1035, an active element which contains three differential 
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amplifiers with emitter follower outputs. Fig. 4.7 gives the 
schematic diagram of the amplifier and logic circuitry. 
Boards containing 16 such units were obtained from the NRC 
Ole themes POL A. —The> logic circuitry had? to be slightly 
modified to suit our needs. Minor changes were also 
introduced to the amplifier circuit to adjust the detection 
threshold above noise level and to reduce crosstalk between 
neighbouring units. 

The first two stages of the active element (a,b) have a 
gain of approximately 72 and are connected as a preamplifier 
to the trigger circuit which makes use of the third stage 
(c). The trigger sensitivity can be adjusted by varying the 
standing current in the diode Ty) by altering Rig To 
eliminate cross talk between neighbouring units, the diode 
current was set at 1.9 mA reducing the overall sensitivity 
of the amplifier to about 1 mV. 

Associated with each amplifier and trigger we have a 
gate Gia) which, in its normal state is low, and a bistable 
memory M(b,c). The additional gate (d) is maintained open in 
our application. Such a gate could be implemented to read 
the state of the memory M by a pulse from the read-out logic 
along with the strobe and reset signal. The trigger output 
is differentiated and inverted by the transistor f5.- The 
pulse crosses gate G changing the state of the memory M. The 
outputs of the 16 bistables are "OR" tied in the encoder, | 


which provides the strobe and reset pulses. 
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The output of the logic circuitry associated with each 
amplifier is fed via an ECL-to-TTL level translator (MC 
10125) into a Fairchild 9318 8-input priority encoder, to 
generate the address of the active input (Fig. 4.86). The 
priority encoder has eight active low inputs (Ig - Ij) with 
highest priority assigned to L. The output (Aga A ) isa 
linearly weighted code representing the active element with 
highest priority. Thus, whenever more than one input is low, 
only the highest priority signal is encoded and the others 
neglected. The inputs are "OR" tied ~OCCTHCIMELO §VYLeLld fa 
group signal (GS) which indicates the presence of an active 
input. An enable input (EI) controls the state of all the 
inputs and overrides them. Whenever FI is high, ail inputs 
are forced high, otherwise they can change state. The enable 
output EO stays high when EI is high or because one of the 
inputs is active low. The proper connection of EO - RI 
allows series or parallel expansion to larger systems. 
Fig. 4.8 shows a series expanded encoder as described in the 
Fairchild semiconductor application handbook’6 and used for 
the 16-wire chamber. 

The time sequence of the signals which controls the MCA 
memory cycle and the outputs of the 16 logic circuits on the 
amplifier board is shown in Fig. 4.9. The sequence begins 
with the -differsntiated trigger output (DTO) at POLL Phun 
Fig. 4.7. The str«cbe and reset signals are generated in the 


encoding circuit following the detection of an event. As 
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FIG. 4.8 Encoder with strobe, reset and start memory signals. 
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FIG. 4.9 Timing sequence of events. 
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Shown in the schematics of Fig. 4.8, the encoder group 
Signal (GS) is ted into a monostable multivibrator (1/2 of 
SN74123). The outiut of one yields the Strobe Signal which 
serves to block additional wires from being encoded until 
the initial event has been stored in the memory. The output 
Q of the second triggers the Start Memory command, while Q 
triggers the Reset Signal. The width of the Strobe Signal 
{ 24 microseconds ) represents the dead time of the system; 
it is due mainly to the timing of the memory cycle. 
Translation from saturated logic to emitter coupled logic 
for both reset and strobe signal is accomplished with a MC 
1017 as shown in Fig. 4.8. 

The encoding system has the apparent disadvantage of a 
non linear response due to the priority scheme of its 
inputs, since given two simultaneous active inputs, only the 
one with the highest priority will be encoded. However, this 
non linearity may be seen to be negligible when the time 
spread bet ween two simultaneous pulses, about 30 


nanoseconds, is compared to the dead time of the systen. 


4.3.3 Data collection 

The encoded signal is stored directly in the memory of 
an ORTEC 6240 MCA using the 5-bit word (GS ApA) 4 Az) - The 
GS bit, which is active for any input, is added to the 
encoded 4-bit word (A oA, A A) to provide 16 odd numbered 


addresses from 1- to 31 ain order that the O channel is 


reserved for timing purposes. The memory cycle is initiated 
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by the Start Memory Signal which is made available along 
with the encoded signal. The cycle itself lasts 20 
“yeaah a read signal counts the level previously 
stored in the particular address, adds one count, then 
writes the revised count back into the memory. After the 
event has been recorded the system is ready to store the 
next event. The data collection capability is better than 
u0-K counts per second. The system was tested by simulating 
the counter pulse using a pulse generator which fed pulses 
into the input of the amplifier board at a repetition rate 
of about 35-K counts per second. The maximum error observed 
in the number of stored events over repeated measurements 
was less than 0.2%. Such rates, are very seldom achieved in 
the detection of labelled proteins movement along the nerve. 
However, the system could be used for other measurements 
such as thin layer chromatography, in which case such 


counting rates could be observed. 


tse! 
















Jas 
' oe : 

pay wok O@2. sere Fea aye, @. tG@b< 
i smushas ele wy «6a AL Bee 
inte: sHee® Mevites wes ne 
ais ES ote 

‘i *o>e to 2908-0 @Ae . Saeee =. 

ad eo ebaee-seq aréees aor Fd - 

’ acon wale @ oebee sabiee. tocapca oad. _ 

: | ret Tddgoe ons Re, Jepas ats fake _ 
1? foo Ge asaue 4 ett opoée 6 
stueo Sesche Fp tna ode Qe eee 








-¢ .~& te: dove. set Oweee eed new 

warts» suet Doping Qi lied. Sipe aot toate wat aa 

st hives omeogaesis TeekeeH oa 

TA, ae ve cate BaD waged nid?: ony! ‘em e 


, oe 84 eee +5 aria augs3 Rr - 
a 
‘ey m 


Wa 


7 seal ’ y <- 


7 








CHAPTER 5 


RESULTS AND DISCUSSION 


To demonstrate the intended applications of the system, 
two types of measurement were performed. In the first, the 
detector was used to record the activity distribution along 
the sciatic nerve of a frog ( Xenopus laevis ); in the 
second type of measurement, it was employed to determine the 
radiochemical purity of a (35S Jj-Methionine source. Fixing 
the gas flow at 80 cc/min and the operating voltage at 1150 
volts, the system was tested for counting efficiency and 
background. A 1*C source was placed against the detector 
window at the center of each of the 16 counters, after which 
the count rate was compared with the known activity of the 
same source. The average value of the efficiency measured in 
such a way waS approximately 70%, with a maximum deviation 
of 5%. Background measurements for the 16 channels were 
slightly higher than expected and ranged between 4 - 8 
cpm/channel. 

In the first type of measurement, the nerve preparation 
(as described section 1.3) was incubated with 275 
microcuries of (35S J-Methionine and 2.5 ml of Ringer's 
solution placed in the spinal cord compartment (refer to 
Fig. 1.2). After an incubation period of approximately 18 


hours, the nerve was severed 15-20 mm distal to the spinal 
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cord and immersed in fresh Ringer's solution in order to 
remove any soluble radioactive material which might have 
leaked from the spinal cord housing into the nerve 
compartment. The nerve was then laid in a tray containing 
Ringer's solution and gently pressed down, by coarse gauze 
stretched across a lightweight plastic frame, against 10 
micrometre thick aluminized mylar. The tray was finally 
positioned on the detector. After the measurements, 
consisting of a set of three 2000-second runs, the nerve was 
removed from the tray and cut in small sections to 
correspond with each counter location. One additional 
measurement with the tray and Ringer's solution alone 
yielded the data for background correction. The nerve 
sections, placed in separate vials, were each dissolved in 1 
ml of Protosol. Following the addition of 10 aml of 
Econofluor, the vials were dark adapted for several hours 
and counted in a Nuclear Chicago Mark I liquid scintillation 
counter. Quenching was uniform among the sections. 

The histograms in Fig. 5.1 illustrate the activity 
distribution along the nerve as measured with the detector 
(Fig. 5.1a), obtained by using the last 30 minute run, and 
the analogous distribution obtained by liquid scintillation 
counting (Fig. 5.1b). The overall shape of the two 
histograms is similar, with the odd displacement in 
configuration from one channel to a neighbouring one. In the 
process of moving the nerve and sectioning it, there is some 


room for error. It is almost impossible to duplicate the 
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(a) 
v8 THO) MLE AB ANSE ALS 15. 16 
COUNTER NO. 
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3 
ae 
3 
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FIG. 5.1 Activity distribution in the nerve obtained with 
a) B -detector, (b) liquid scintillation counter. 
The total activity in (a) amounts to0.6% of (b). 
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relative position of the segments with respect to the 
multiple detector when the nerve is removed and segmented. 
This means that in some cases, liquid scintillation counting 
is being performed on a segment whose length has been 
changed (hence change in activity) or on one that does not 
correspond identically to its position in the detector. 
Therefore, the position of the nerve seen by any given 
channel may not be the same as seen in liquid scintillation 
counting. This would account for the displacement in the 
major peak of the two histograms in Fig. 5.1. 

In the second type of measurement performed with the 
system, we tested the radiochemical purity of the S-emitter 
used in the nerve experiment, (35S J-Methionine. A 
chromatographic separation was used which allows 
identification and quantification of the radiochemical 
inpurities present in the labelled amino acid. The principle 
employed in chromatography is based on the different 
migration velocities of the components of the sample in the 
mobile phase which is provided by the solvent. The 
stationary phase consists of an absorbent. Measurements were 
performed using a 10-cm strip of ITLC-SG chromatographic 
paper in a Seprachron chamber both of which were 
manufactured by Gelman Instrument Company, together with a 
solvent whose composition of n-Butanol:NH,OH:H,0 (30:1: 2) 
achieved optimum separation of the radiochemicals. The 
chromatographic strip was first spotted with 1 microlitre of 


(35S ]-Methionine and then placed inside the Seprachrom 
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chamber containing the solvent. As the solvent rises through 
the medium by absorption and capillary action, the various 
chemical components will rise with different velocities. 
Methods of assaying the radioactivity distribution in the 
paper strip makes use of visual observation employing 
chromogenic compounds and of commercial radiochromatogranm 
scanners. 

Our system offers an alternate method for determining 
the activity distribution, which is achieved by placing the 
chromatographic strip against the detector window for any 
length of time. Fig. 5.2 shows a distribution SMA CLeV LLY 
measured with the present detector. Point A in Fig. 5.2 
denotes the origin, while D indicates the location of the 
solvent front. The peaks at B and C represent the impurity 
of Methionine sulfoxide and the Methionine respectively. The 
percentage of radiochemical purity (%RP) can be determined 
from the number of counts in peak B and peak C from the 


simple relation 


j m counts peak C 
ae UES counts peak C + counts peak B 


In the example shown, the radiochemical purity is 74%. 


ae <a ae Ser Ee ED SESS AE OT Om RD 


The detection system described in chapter 4 was devised 
mainly for the study of axonal transport. It has been shown 


that the system can also be used in chromatographic 
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determination of the radiochemical purity of labelled 
compounds. The detector can be easily assembled and allows 
for quick replacement of any part. The read-out system, 
adopting the one amplifier per wire and priority encoder, 
offers a simple method of extracting the information. 
Because of the fact that the system is still in its 
development stage, there are still certain inherent 
difficulties that have to be dealt with, before the 
construction of a larger unit. Most significant of these is 
crosstalk, which originates in the amplifier and trigger 
circuit and in the connector used between the counter and 
the high voltage network. In the latter case, a connector 
with smaller capacitance between pins would suffice to 
eliminate such crosstalk. More troublesome, however, is that 
occurring between neighbouring units in the amplifier, 
trigger and memory board. Indications are that this 
crosstalk occurs at the trigger stage because of poor 
grounding between units. With the present system this source 
of crosstalk was greatly reduced by lowering the sensitivity 
of the trigger circuit; however, this was accompanied by a 
loss in detection efficiency. The crosstalk level between 
the 16 amplifier units, although reduced, remains at about 


4-8%. 
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APPENDIX A 


———— SS ee ee ee SaaS Se Se 


The proportional counter of rectangular geometry 
consists simply of a central electrode of circular cross 
section enclosed in a long hollow rectangular metal box 
usually kept at ground potential. The counter cross section 
is assumed to be small compared to the counter length so 
that the potential satisfies Laplace's equations in two 
dimension. 

The analytical solution has been given by Tomitani®2 
using the method of conformal transformation by means of 
elliptic sine functions. With such a technique one obtains 
the solution to the electrostatic problem both for a single 
anode and for a series of anodes placed periodically within 
the cathode, and is applicable to those problems which have 
periodicity in two directions. 

For anode-cathode arrangements which are not 
symmetrical, a numerical approach is advantageous as 
compared to this type of analytical approach. This permits 
application to more complex anode-cathode configurations. 
One such numerical method makes use of Fredholm's first 


integral equation to solve the boundary value problem®*, 


The axial symmetry of the counter reduces the 


electrostatic problems to a two-dimensional one. Assume the 
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following: 

L is the outer boundary delimited by the counter, 

R is the radius of the anode wire, 

D is the domain bounded externally by the counter 
(cathode) and internally by the anode, and 

V is the potential on the anode. 

In addition, it is assumed that R is small compared to 
the counter dimensions : In such a counter, the 
electrostatic potential, ¢ , in D satisfies the 


two-dimensional Laplace equation 





with the boundary conditions 


0 on L 


(x,y) = 


V on the anode surface 


Equation A.1 may be rewritten in integral form as¢S 
> > > > aa > 
o@) = f o(q) log \a- >| dq + A log |B - w| , A.2 
L 


where the integral is taken over the boundary L, Pp is the 
position vector of a point within D, q is the position 
vector of a point on L,and w is the position vector of the 


centre of the anode. The problem is solved by determining 


Aq 
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old), the charge density on L , and the constant A. Imposing 
the boundary condition on L, i.e. 


o(p) = 0 


equation A.2 yields 
> > 
Jf 0@) 108 \a- >| dq +A log |p -w| = 0 ee 
L 


for p on the boundary L. The second boundary condition, on 


the anode surface, becomes 


> > > 
J 0G we |3-%] aq +a toe R=v,=1 Det 
ile 


a 


The replacement V, = 1 is to simplify the calculations. 
Equations A.3 and A.4 have to be solved simultaneously for 
O(q) and A. 

Consider now a rectangular counter as in Fig A.1. The 
boundary L is defined by the path OABCO; L is divided into n 
segments, each having a charge distribution Oj} (j=1,---,N). 
If the segments are small enough, it may be assumed that 
O. is constant on each segment. The integral over L (in 


J 


equation A.3) can then be replaced by a sum of terms, 
e > > > > 
Sy 0, if log | a a, | dd teas Loe | ah, = ly | = 0 A.5 
jai j ak i 


for i=1,.+e,ne For n intervals (or segments), there will be 
n linear equations. 

On the anode surface (equation A.4), a similar approach 
yields the following 


uw > > 
ie 0, f roe [a-% | aq + arog R= 1 A.6 
Uae j 
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FIG. A.1 Rectangular counter with anode at (a, 8). 
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Equations A.5 and A.6 form a set of (n + 1) simultaneous 
equations in the (n + 1) unknowns oj and A. The solution to 
the Laplace equation then becomes 


~ > 
GAP 


n 
Pai e ao al log 


j 














where n indicates the number of segments. The approximation 
of % (p) (Equation A.2) by op) becomes more reliable with 


increasing n. 


A.3 Results 

The potential ¢ (p) is obtained by first solving 
equations (A.5) and (A.6) Simultaneously. The solutions OF 
and A are then substituted into equation (A.7) to yield 
op). This whole procedure is done utilizing the computer 
program MAT. The solutions 6” (p) are then plotted as 
equipotentials, varying from 5% to 60% of anode voltage, in 
IIcCrements (Of o% ~(Or,, with aVae = 1, from 0.05 to, 0.6 in 
increments of 0.05). The computing was done on the 
University of Alkerta Amdahl 470-V6 computing facility. The 
numerical values for > (Pp) were contour plotted using the 
Calcompq plotter. A listing of the entire computer program 
follows in the next section. 

Results from two counter configurations were obtained. 
These particular configurations were chosen so as to enable 
comparison with the results of Tomitani. The eaioocent aie 


obtained in this manner, for a Square counter, are shown in 
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FIG. A.2 Square counter. Equiportential lines 
obtained a) with numerical technique 
b) with exact solution up to 30% of anode 


voltage. 
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(b) 


FIG. A.3 Rectangular counter. Equipotential lines obtained a) with 
numerical technique b) with exact solution up to 30% of 
anode voltage. 
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With analogous results using the exact 


similar comparison for a rectaagular 


counter is demonstrated in Fig. A.3. These results show that 


the differences between 


techniques, 


A.4 Program "MAT" 


200 
100 


equipotential lines, in the two 


become smaller near the cathode. 


= 


DIMENSION A1(20,20),A2(20,20) ,A3 (20,20) ,B1 (20, 20) 


DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 


DD1¢21,21) ,DD2(21,21) ,DD3 (21,21) ,ADD1 (21, 21) 
ADD2 (21,21) ,ADD3(21,21) ,ADD4 (21,21) 
AD1(21,21) ,AD2 (21,21) ,AD3 (21,21) 

AD4 (21,21) ,D1(21,21) ,D2(21,21) ,D3(21,21) 

D4 (21) ,D5 (21) ,DD5(21) 

C1 (20,20) ,C2 (20,20) ,C3 (20,20) , AAT (20, 20) 
AA2(20,20) ,BB2 (20,20) ,CC2 (20,20) , 8A 3( 20,20) 
AB(81,81) ,B2 (20,20) ,B3 (20,20) ,BB1(20, 20) 
CC1(20,20) ,BB3 (20,20) ,CC3 (20,20) 

A4 (20) ,B4 (20) ,C4 (20) ,A5 (20) , B5 (20) ,C5 (20) 
AA5 (20) ,BB5(20) ,CC5 (20) 

VR(65) ,MA (65,65) 

CM (50,50) ,CVAL(12) , LOP (7) , VOP (7) 


DATA CVAL/12*0.0/, IOP/7*0/,CM/2500*0.0/ 

REAL L,LL,MA 

WRITE(6,10) 

FORMAT(*ENTER VALUES OF L,LL,RX,RY(=LOCATION OF ANODE) * 
READ(5,20) L,LL,RX, RY 


FORMAT (F6. 


4) 


ANODE RADIUS 


R=0.01016 


WRITE(6,30) 


FORMAT(*ENTER VALUES OF N,JR, 


N=16 JR=1") 


READ(5,40) N,JR 
FORMAT (213) 


M=J R*N 
NN=N+1 


DO 100 I=1,NN 
DO 200 K=1,N 


J=I-1 


A=FLOAT (J) *L/FLOAT (N) -FLOAT (2*K-1) *L/FLOAT (2*N) 
A1(K,1) =A*ALOG (ABS(A)) -A 


CONTINUE 
MM=M+1 


DO 110 I=1,M 
DO 210 K=1,N 


J=I-1 
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1106 


220 
120 


1007 
1006 


400 
300 


410 
310 


420 
320 


1009 
1008 
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B=L-FLOAT (2*K-1) *L/FLOAT (2*N) 
C=PLOAT (J) *LL/FLOAT (M) 

B1(K,1)=C/2*ALOG (B¥*2+C**2) -C+B¥ATAN (C/B) 
CONTINUE 

NN=N41 

DO 120 I2=1,NN 

I=NN+1-12 

DO 220 K=1,N 

J=I-1 

D=FLOAT (J) *L/FLOAT (N) ~ FLOAT (2*K-1) *L/FLOAT (2*N) 
C1(K,1) =D/2.0*ALOG (D** 2+LL**2) -D+LL*ATAN (D/LL) 
CONTINUE 


MM=M+1 

DO 1006 I2=1,MM 
T=AM+ 1-22 

DO 1007 K=1,N 
J=I-1 


A=FLOAT (J) *LL/FLOAT (M) 
B=FLOAT (2*K-1) *L/FLOAT (2*N) 

D1 (K,1) =A/2*ALOG (B*¥ *2+A**2) -A+B¥*ATAN (A/B) 
CONTINUE 

NN=N+1 

DO 300 I=1,NN 

DO 400 K=1,M 

J=I-1 

A=FLOAT(J) *L/FLOAT(N) -L 

B=FLOAT (2*K-1) *LL/FLOAT (2*M) 

A2 (K,1) =A/2.0*ALOG (A**2+B**2) - A+ B¥ATAN (A/B) 
CONTINUE 

MM=M+1 

DO 310 I=1,MM 

DO 410 K=1,M 

J=I-1 

A=FLOAT (J) *LL/FLOAT (M) -FLOAT (2*K-1) *LL/FLOAT (2*M) 
B2 (K, 1) =A*ALOG (ABS (A)) -A 

CONTINUE 

NN=N+1 

DO 320 I2=1,NN 

I=NN4+1-12 

DO 420 K=1,M 

J=I-1 

A=FLOAT (J) *L/FLOAT (N) -L 
B=LL-FLOAT (2¥*K-1) *LL/ FLOAT (2*M) 

C2 (K, 1) =A/2.0*ALOG (A*¥*2+B**2) ~A+B*ATAN (A/B) 
CONTINUE 


MM=M+1 

DO 1008 I2=1,MM 
I=MM+1-I2 

DO 1009 K=1,M 
J=I-1 


A=FLOAT (J) *LL/FLOAT (M) -FLOAT (2*K-1) *LL/FLOAT (2*M) 
D2 (K,1I) =A/2*ALOG (A*¥*2+L**2) -A+L*ATAN (A/L) 
CONTINUE 
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600 
500 


610 
510 


620 
520 


10:41 
1010 


1013 
1012 


NN=N+1 

DO 500 I=1,NN 

DO 600 K2=1,N 

K=N+1-K2 

J=I-1 

D=FLOAT (J) *L/FLOAT (N) - FLOAT (2*K-1) *L/FLOAT (2*N) 
A3(K,1) =D/2.0*ALOG (D** 2+LL**2) -D+ LL*¥ATAN (D/LL) 
CONTINUE 

MM=M+14 

DO 510 I=1,MM 

DO 610 K2=1,N 

K=N+ 1-K2 

J=I-1 

A=FLOAT (J) *LL/FLOAT (M4) -LL 
B=L-FLOAT (2*K-1) *L/FLOAT (2*N) 

B3 (K, 1) =A/2.0*ALOG (A**2+B**2) -A+B*ATAN (A/B) 
CONTINUE 

NN=N+1 

DO 520 I2=1,NN 

I=NN4¢17-1I2 

DO 620 K2=1,N 

K=N+t+ 1-K2 

J=I-1 

A=FLOAT (J) *L/FLOAT (N) -FLOAT (2*K-1) *L/FLOAT (2*N) 
C3(K,1) =A*ALOG (ABS (A))-A 

CONTINUE 

MM=M+1 

DO 1010 I2=1,MM 

I=MM41-1I2 

DO 1011 K2=1,N 

K=N+ 1-K2 

J=I-1 

A=FLOAT (J) *LL/FLOAT (M) -LL 

B=FLOAT (2*K-1) *L/FLOAT (2*N) 

D3(K,1) =A/2*ALOG (A**2+B**2) -A+B¥*ATAN (A/B) 
CONTINUE 

NN=N+14 

DO 1012 I=1,NN 

DO 1013 K2=1,M 

J=I-1 

K=M+ 1-K2 

A=FLOAT (J) *L/FLOAT (N) 

B=FLOAT (2*K-1) *LL/FLOAT (2*M) 

AD1(K, 1) =A/2*ALOG (A**2+B**2) -A+B*ATAN (A/B) 
CONTINUE 

MM=M+1 

po 1014 I=1,MM 

DO? 1015) Kz=1, 

J=I-1 

K=M+ 1-K2 

A=FLOAT (J) * LL/FLOAT (M) 

B=FLOAT (2*K-1) *LL/ FLOAT (2*M) 
ZEN=AT AN ( (A-B) /L) 


108 





oe, nee areas 
f+ 


roan ory S sa ZrL) Casey a 





















~, 


eee gy CP 


- fet O0¢.08 
doce (19.08% 
s%-F hag. 9S? 
h-yah . 
hp te Oot heeg tg T4029 88 
ae  £~Fee tA - ee a 
ad h ©1200 44 4° a. 388 Gre; - = 
aor Taos, C76 


rohe@s i > : 
24 sy 90% Oa aon 
-f Bet 0 : 
dbo “a we wae a“ 
2 f ‘tod - : nar 
. ori’ : Our 
week a 8 ae | 
(~ (07 oe Sage eee? oh a 
Telteo DSR 
iehehe i” 
r,*o6 Cl ober Ge Pa ’ 


‘ le rehnet "I f 


sitege “c8r, oe, ,i4 
Sie pkey ia 7 
i =i eh ' j 


af ti y's 495 tae bp Open’ 


ave i ua 


iT grel t 
o.Fae, tf 





7 
Pine 


» = 


109 


1015 AD2 (K,1)=(A-B) /2*ALOG (L**2+ (A-B) **2) - (A-B) +L*ZEN 
1014 CONTINUE 

NN=N+1 

DO 1016 I2=1,NN 

I=NNt+1-I2 

DO 1017 K2=1,M 

K=M+1-K2 

J=I-1 

A=FLOAT (J) *L/FLOAT (N) 

B=LL-FLOAT (2*K-1) *LL/FLOAT (2*M) 
1017 AD3(K, I) =A/2*ALOG (A**2+B**2) -A+B*ATAN (A/B) 
1016 CONTINUE 

MM=M+1 

DO 1018 IT2=1,MM 

I=MM41-12 

DO 1019 K2=1,M 

K=M+ 1-K2 

J=I-1 

A=FLOAT (J) * LL/FLOAT (M) 

B=FLOAT(2*K-1) *LL/ FLOAT (2*M) 
1019 AD4 (K, 1) = (A-B) *ALOG (ABS (A-B) ) - (A-B) 
1018 CONTINUE 

DO 201 K=1,N 

DO 101 IT=1,N 

J=I+1 

AA1(K, 1) =A1(K,J)-A1(K, I) 
101 CONTINUE 
20 1 CONTINUE 

DO 811 I=1,N 

DO 812 J=1,N 
812 CONTINUE 
811 CONTINUE 

DO 211 K=1,N 

DO 111 I=1,8 

J=I+1 
171 BB1(K, 1) =B1 (K, J) -B1 (K, 1) 
211 CONTINUE 

DO 221 K=1,N 

po 121 I=1,N 

J=Nt2-TI 

JJ=Nte1-1 
121 CC1(K, I) =C1(K, J) -C1(K,JJ) 
ee CONTINUE 

DO 1020 K=1,N 

DO 1021 I=1,M 

J=M+2-I 

JJ=Mt+1-1 
1021 DD1(K,1)=D1 (K,J)-D1(K,JJ) 
1020 CONTINUE 

DO 401 K=1,M 

DO 301 I=1,N 

J=I+t+1 

AA2 (K, 1) =A2 (K, J) ~A2 (K, 1) 
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CONTINUE 

CONTINUE 

DO 411 K=1,M 

DO 311 I=1,M 

J=I+1 

BB2 (K, 1) =B2 (K,J) -B2 (K, 1) 
CONTINUE 

DO 421 K=1,M 

DO 321 I=1,N 

J=N+2-1 

JJ=N41-T 

CC2 (K, I) =C2 (K, J) -C2 (K, JJ) 
CONTINUE 

DO 1022 K=1,M 

DO 1023 I=1,M 

J=M+2-1 

JJ=M+1-1 

DD2 (K, I) =D2 (K, J) -D2 (K, JJ) 
CONTINUE 

DO 601 K=1,N 

KK=N+1-K 

DO 501 I=1,N 

J=I+1 

AA3 (K, I) =A3 (KK,J) -A3 (KK,1) 
CONTINUE 

DO 611 K=1,N 

KK=N+1-K 

DO 511 I=1,M 

J=I+1 

BB3 (K, I) =B3 (KK,J) -B3 (KK, 1) 
CONTINUE 

DO 621 K=1,N 

KK=N+1-K 

DO 521 I=1,N 

J=N+2-I 

JJ=N+1-I 

CC3 (K, I) =C3 (KK,J) -C3(KK,JJ) 
CONTINUE 

DO 1024 K=1,N 

J=N+1-K 

DO 1025 I=1,M 

JJ=M+2-I 

JIJ=M+ 1-1 

DD3 (K, 1) =D3 (J, JJ) -D3 (J, JJJ) 
CONTINUE 

DO 1026 K=1,M 

J=M+1-K 

po 1027 I=1,N 

JJ=I+1 

ADD1 (K,1I) =AD1(J,JJ) -AD1(J,1I) 
CONTINUE 

DO 1028 K=1,M 

J=M+1-K 





DO 1029 I=1,M 
ies 
1029 ADD2 (K,1I)=AD2 (J, JJ) -AD2(J,1) 
1028 CONTINUE 
DO 1030 K=1,M 
J=M+1-K 
DO 1031 I=1,N 
JJ=N+2-1 
JJJ=N+1-1 
1031 ADD3 (K,I) =AD3(J,JJ) -AD3(J,JJJ) 
1030 CONTINUE 
DO 1032 K=1,M 
J=M+1-K 
DO 1033 I=1,M 
JJ=M42-1 
JIT=M+1-1 
1033 ADD4 (K,1I)=AD4(J,JJ) -AD4 (J, JJJ) 
1032 CONTINUE 
DO 250 I=1,N 
DO 251 J=1,N 
251 AB(I,J}=AA1 (I,J) 
DO 252 J=1,4 
252 AB(I,J+N)=BB1(I,J) 
DO 253 J=1,N 
253. AB(I,J+N#M)=CC1(I,J) 
DO 1035 J=1,M 
1035 AB(I,d+Nt¢M+N)=DD1 (I,J) 
250 CONTINUE 
DO 260 I=1,M 
DO 261 J=1,N 
261 AB (N+I,J)=AA2 (I,J) 
DO 262 J=1,M 
262 AB(N+I,J+N)=BB2(I,J) 
DO 263 J=1,N 
263 =AB(NtI, J+N+#M)=CC2 (I,J) 
DO 1034 J=1,M 
1034 AB (N+I,J¢N¢M+N)=DD2 (I,J) 
260 CONTINUE 
DO 270 I=1,N 
DOe27 jed=4 7h 
271 AB(N+M41I,J)=AA3 (I,J) 
DO 272 J=1,¥ 
272. = AB(N#M+I,J+#N) =BB3 (I,J) 
DO 273 J=1,N 
273 AB(N4#M+I, J+N¢M) =CC3 (I,J) 
DO 1036 J=1,M 
1036 AB(NtM41,J+N+M+#N) =DD3 (I,J) 
270 CONTINUE 
DO 1037 I=1,M 
DO 1038 J=1,N 
1038 AB(N+M4+N4+I,J)=ADD1 (I,J) 
DO 1039 J=1,M 
1039 AB(N¢M+N+I,N+J)=ADD2 (I,J) 
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1040 


1041 
1037 


Tig 


888 


999 


1001 


778 


7179 


780 


1002 


889 


890 


89 1 


DO 1040 J=1,N 

AB (N4#M#N¢1,N¢M+J) =ADD3 (I,J) 
DO21041 /J=1,é 

AB (N*¢M#N4I, N¢M#N¢J) =ADD4 (I,J) 
CONTINUE 

AL=RX*L 

ALL=RY*LL 

DO 777 I=1,N 
A=FLOAT (2*1I-1) *L/FLOAT (2*N) 
A4(T)=0.5*ALOG ((A-AL) **24ALL**2) 
DO 888 I=1,M 

A=FLOAT(2*1I-1) *LL/FLOAT (2*M) -ALL 
B4 (1) =0.5*ALOG ((L-AL) **2+A**2) 
DO 999 I=1,N 

J=Nt+ 1-1 

A=FLOAT(2*J-1) *L/FLOAT (2*N) 

C4 (1) =0.5*ALOG ((A-AL) **2+4 (LL-ALL) **2) 
DO 1001 I=1,M 

J=M+1-T 

A=FLOAT (2*J-1) *LL/ FLOAT (2*M) ~ALL 
D4 (1) =0.5*ALOG (AL**2+A **2) 
NN=N+1 

DO 778 I=1,NN 

J=I-1 

A=FLOAT (J) *L/FLOAT (N) -AL 

A5 (I) =(A/2.0) *ALOG (A¥*2+ALL**2) -A+ALL*ATAN (A/ALL) 
CONTINUE 

MM=M+1 

DO 779 T=1,MM 

J=I-1 

A=FLOAT (J) *LL/FLOAT (M) -ALL 

B5 (1) =A/2.0*ALOG( (L-AL) **2+A**2) -A+ (L-AL) *ATAN (A/ (L-AL) 
NN=N+1 

DO 780 I=1,NN 

J=NN-I 

A=FLOAT (J) *L/FLOAT(N) -AL 
ZEN1=ATAN=(A/ (LL-ALL) ) 

C5 (1) =A/2.0*ALOG (A**2+ (LL-ALL) **2) ~A+ (LL~ALL) *ZEN1 
MM=M+14 

DO 1002 I=1,MM 

J=MM-I 

A=FLOAT (J) *LL/FLOAT (M) -ALL 

D5 (1) =A/2*ALOG (AL**2+A**2) -AtAL*ATAN (A/AL) 
DO 889 I=1,N 

J=I+1 

AAS (1) =A5(J)-A5 (I) 

DO 890 I=1,mM 

J=I+1 

BB5 (I) =B5 (J) -B5 (I) 

DO 891 I=1,N 

J=I+1 

CC5 (1) =C5(1)-C5 (J) 

DO 1003 I=1,M 
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J=I+1 

DD5 (1) =D5 (I) -D5(J) 
N1I=N¢M+N+M 

DO892 J=1,N 

AB(N1+1,J) =AA5 (J) 

DO 893 J=1,M 
AB(N1+1,N+J) =BB5 (J) 

DO 894 J=1,N 
AB(N1+1,N+M+J) =CC5 (J) 
DO 1004 J=1,4 

AB(N1+1, NeM4*N+J) =DD5(J) 

DO 895 I=1,N 

AB(I,N1+1) =A4 (1) 

DO 896 I=1,M 
AB(Ntl,N1+1)=B4(T) 

DO897 I=1,N 

AB(NtM41I,N14+1) =C4 (1) 

DO 1005 T=1,M 

AB( N#MtN+1,N14+1) =D4 (1) 
N2=N1+1 

AB(N2,N2) =ALOG (Rk) 


COMPUTATION OF THE VALUES OF SIGMA AND A 


WRITE(6, 1) 
FORMAT(* ENTER VALUES FOR GELG IM,IN,IESP') 
IM=2¥*M+2*N+1 
IN=1 
TESP=1000000 
DO 3 J=1,N2 
DO 4 IT=1,N2 
MA (I,J) =AB(I,J) 
CONTINUE 
CONTINUE 
WRITE(6,5) 
FORMAT('ENTER VALUES OF VECTOR VR‘) 
N3=N2-1 
DO 9 I=1,N3 
VRq«1)=0.0 
VR(N2) =1.0 
CALL GELG(VR,MA,IM, IN, LESP, IER) 
WRITE(6,7) 
FORMAT (*SOLUTION OF LINEAR EQTS.'°) 
WRITE(6,8) (VR(I) ,L=1,1M) 
FORMAT (' ¢,1X,8E15. 7) 
WRITE (6,25) 
FORMAT(*ENTER NO. OF GRID LINES IN X AND Y 
READ (5,26) I1I,JJ 
FORMAT (213) 
DXx=0.0 
DY=0.0 
DO 444 K1=1,dJI 
DY=LL/FLOAT (JJ) *F LOAT (K1) 


II,JJ,11t=JJ 





DO 445 K2=1,II 
E1=0 
DX=L/FLOAT (II) *FLOAT (K2) 
DO 333 J=1,N 
A=FLOAT (J) *L/FLOAT (N) -DX 
B=FLOAT (J-1) *L/FLOAT (N) -DX 
IF(DY.EQ.0.0)GO TO 448 
C=A/2*ALOG (A¥* 2+ DY¥**2) -A-DY*ATAN (A/ (-DY) ) 
D=B/2* ALOG (B¥* 2+DY**2) -B-DY*ATAN (B/(-DY)) 
GO TO 409 
448  IF(A.EQ.0.0) GO TO 447 
C=A/2*ALOG (A**2) -A 
446 IF(B.eEQ.0.0) GO TO 449 
D=B/2* ALOG (B**2) -B 
GO TO 409 
447 c=0.0 
GO TO 446 
449 D=0.0 
409 E=VR(J)*(C-D) 
333 EI=E1+E 
E2=0 
DO 334 J=1,M 
A=FLOAT (J) *LL/FLOAT (M) -DY 
B=FLOAT (J-1) *LL/FLOAT (M) -DY 
IF((L-DX).EQ.0.0) GO TO 404 
C=A/2* ALOG ( (L-DX) **2+A**2) -A+ (L-DX) ¥ATAN (A/ (L-DX) ) 
D=B/2* ALOG ( (L- DX) **2+B**2) -B+(L-DX) *ATAN (B/ (L-DX) ) 
GO TO 408 
404 IF (A.EQ.0.0) GO TO 405 
C=A/2*ALOG (A*¥*2) -A 
406 IF(B.EQ.0.0) GO TO 407 
D=B/2*ALOG (B**2) -B 
GO TO 408 
G05 C=0u0 


uO?) b=0.0 
408  E=VR(NtJ)*(C-D) 
334 E2=E2+E 
E3=0 
DO 335 J=1,N 
A=FLOAT (N+1-J) *L/ FLOAT (N) -DX 
B=FLOAT (N-J) *L/FLOAT (N) -DX 
IF((LL-DY).EFQ.0.0) GO TO 503 
C=A/2* ALOG (A** 2+ (LL-DY) **2) -A+ (LL-DY) *ATAN (A/ (LL-DY) ) 
D=B/2*ALOG (B¥* 2+ (LL-DY) **2) -B+ (LL-DY) *ATAN (B/ (LL-DY) ) 
GO TO 509 
503 IF(A.EQ.0.0) GO TO 504 
C=A/2* ALOG (A**2) -A 
505 IF(B.EQ.0.0) GO TO 506 
D=B/2*ALOG (B*¥*2) -B 
GO TO 509 
504 c=0.0 
GO TO 505 
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D=0.0 
E=VR (N+M+J) * (C-D) 
E3=E3+E 
E4=0.0 
DO 1055 J=1,M 
A=FLOAT (M+1-J) *LL/ FLOAT (M) -DY 
B=FLOAT (M-J) *LL/FLOAT (M) -DY 
IF(DX.EQ.0.0) GO TO 1046 
C=A/2* ALOG (DX*¥*2+A**2) -A-DX*ATAN (A/(- DX) ) 
D=B/2*ALOG (DX**2+B**2) -B-DX*ATAN (B/ (-DX)) 
GO TO 1045 
IF(A.EQ.0.0)GO TO 1047 
C=A/2* ALOG (A**2) -A 
IF(B.EQ.0.0) GO TO 1049 
D=B/2*ATAN (B**2) -B 
GO TO 1045 
c=0.0 
GO TO 1048 
D=0.0 
E=VR (N*tM+NtJ) * (C-D) 
E4=E4+t+E 
EO 1=E1+E2+E3+Eu4 
CM(K2,K1) =EO1*VR(N2) /2*A LOG ( (DX-AL) **2+ (DY- ALL) **2) 
CONTINUE 
CONTINUE 
WRITE(8, 2000) ((CM (I,J) ,1=1,K2) ,J=1,K1) 
FORMAT (E15.7) 
CALL PLOTS 
WRITE (6,646) 
FORMAT(*ENTER VALUES OF XSIZE,YSIZE =L,LL') 
READ (5,647) XSTZE,YSIZE 
FORMAT (F4, 2) 
NC=12 
TOP (1) =1 
TOP (2) =1 
IOP (3) =0 
TOP (4) =0 
IOP (5) =0 
TOP (6) =0 
IOP (7) =1 
VALUES OF CONTOUR LINES TO BE PLOTTED 
CSET=0.0 
DO, 27. L=1,12 
CVAL (I) =CSET+0.05 
CSET=CVAL (TI) 
CONTINUE 
T1I=II 
JiJ=JJI 
CALL CONTUR(XSIZE,YSIZE,CM,50,111,J1J3,CVAL,NC,IOP, VOP) 
CALL PLOT(0.0,0.0,999) 
STOP 
END 
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